(§) RfePUBLIQUE FRANCAISE 

INSinUT NATIONAL 
DE LA PR0PR1ETE INDUSTRIELLE 



(fi) N° da publication : 2 591 884 

(ft ninabcr qua pour tes 
conunandes do reproduction) 

<2y N° d'enregistroment national : 86 17673 



PARIS 



(g) Int CI* : A61 B 8/14; G06 F 15/42. 



DEMANDE DE BREVET D'INVENTION 



A1 



Data da depot : 17 decembre 198a 
Priorita : US. 17 decembre 1985. n» 810,872. 



(43) Date da la mise a disposition du public da la 
damanda : BOPI « Brevets » n° 26 du 26 juin 1987. 

(60) References a d'autres documents nationaux appa- 



@ Demandeurts) : WASHINGTON RESEARCH FOUNDA*' 
WK Socm constitute sous les his de rEtat de Was- 
hmgton. ~ U& 



@ Inventeurte) : Francis MamveJ! Greene Jr et Donald Eu- 
gene Strandness Jr. 



ITtulairets) : 



74) Mandataire(s) : Cabinet William J. Rezac. 



Precede non envahissant da diagnostic da la stenose vasculaire. 



(5?) Precede et dispositif permettent de diagnostiquer auto* 
matiquement et objective ment le degre de stenose cfun vais- 
seau a partir de donnees Doppter cornbinees en mode B et a 
impulsions, anaJysees en segments spectra ux dans un anaJy- 
seur de spectres a transformation rapide de Fourier en temps 
reel 

Le dispositif comprend un ordinateur numerique convenable- 
ment programme pour effectuer la reconnaissance de formes 
des segments spectraux ou d'une moyenne (Tensembte des 

< segments spectraux pour classer le degre de stenose en 
fonction d'une base da donnees de stenoses connues en 
I provenance de segments spectraux classes a favance. 

La figure 2 est une representation schematique (fun signal 
^ Doppter type sur une courbe {Tampfitude en fonction du 
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PROCEDE NON. ENVAHISSANT DE DIAGNOSTIC DE LA 
STENOSE VASCULAIRE 
Cette invention concerne ion proc^de automatique, objec- 
tif et non envahissant d 1 interpretation du degr£ de st£nose 

5 (r^trecissementjd'un vaisseau & partir des donnees de varia- 
tion de frequence (due h I'effet Doppler) ultrasonore de 
la circulation du sang dans le vaisseau. Plus parti culiferement, 
1' invention concerne un proc6d6 empirique inforraatique de 
reconnaissance statistique des formes entre les donn^es du 

10 patient et une base de donnees de spectres connus et classes 
a l'avance, pour diagnostiquer le degr£ de st£nose des art&res 
carotides d'un patient en di verses categories de st^nose. 
Le procede fait appel & une scanographe a ultrasons, un ana- 
lyseur de spectre a transformation rapide de Fourier , et 

15 un ordinateur num^rique specialise pour ef f ectuer . 1 'analyse 
des donnees provenant des arteres carotides primitive et 
interne. 

La detection et !• Evaluation de 1 'arteriopathie caroti- 
dienne sont importantes pour le traitement des ischdmies 

20 cer6brales transitoires et pour la prevention des attaques* 
L'art^riographie, methode envahissante et cofiteuse, pr^sente 
des risques de morbidity et de mortalite qui limitent son 
emploi h des patients s£lectionn£s . On met actuellement au 
point un procede non envahissant de diagnostic h. ultrasons 

25 en tant que solution de rechange & l'art^riographie. 

Pour la pathologie d f occlusion des artferes carotides 
qui est associ£e aux anomalies de l'ecoulement sanguin, telles 
qu'elles sont detect^es avec un appareil a ultrasons & effet 
Doppler, des diagrammes d'dcoulement irreguliers proches 

30 de la zone l£see ont 6t6 6tud±6s par H. Keller, et al . , "Non- 
invasive Measurement of Velocity Profile? and Blood Plow 
in the Common Carotid Artery by Pulsed Doppler Ultrasound," 
Stroke , Juil.-AoQt 1976 r 7(4), 370-7, et ont 6t6 en outre 
quantifies par R. Rutherford et al. , "The Use of Velocity 

35 Waveform Analysis in the Diagnosis of Carotid Artery Occlusive 
Disease," Surgery 695-702 (1977), qui utilisaient des para- 
mfetres de forme d'onde mesures & la main. Des parametres 
mesures par ordinateur, obtenus a partir de l'art&re carotide 
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primitive, ont 6te utilises par F. Greene, Jr., "A Micro- 
processor Based Pattern Recognition Approach to Diagnosing 
Atherosclerosis, H Master's Thesis, University of Washington 
(1979) et par S. Sheriff, "Mathematical Feature Extraction 
5 Applied to the Entire Doppler Shifted Frequency Signal Obtained 
from the Common Carotid Artery (Abstract), Blood Flow Theory 
and Practice , The Biol, Eng'g Soc, Londres (1980) 34-1, 
34-4. 

Pour determiner l'etendue de l'art£riopathie h partir 
10 de donn£es de variation de frequence due a I'effet Doppler 
par ultrasons a partir des art&res carotides primitive et 
interne, on utilise le procede assiste par ordinateur de 
la presente invention pour reconnaitre des formes ou des 
structures dans les caracteristiques spectrales et tempe- 
rs relies du signal Doppler retrodif f use. Cette m£thodologie, 
appelee "reconnaissance statistique des formes", a trouve 
application dans des disciplines qui vont de la recherche 
chimique (B. Kowalski, "Pattern Recognition in Chemical 
Research," Computers in Chem. and Biochem. Res. , vol. 2, 
20 Academic Press, N.Y . (1976) a 1* identification d'anomalies 
du foie a l'aide du diagnostic par ultrasons (R. Lerski et 
al. , "Computer Analysis of Ultrasonic Signals in Diffuse 
Liver Disease," Ultrasound Med. Biol. , 5, 341-350 (1979). 
Une telle approche convient trfes bien au problfeme de l'6va- 
25 luation de la st£nose vasculaire, dans les cas oil, tandis 
que la nature precise des processus physiques mis en jeu 
n'est pas bien 6lucidee, on peut encore determiner l'etendue 
du mal. Les algorithmes mis en jeu dans le proc£d£ de la 
presente invention ont permis d'obtenir des r£sultats utiles 
30 du point de vue du diagnostic, avec des donn£es cliniques 
obtenues dans des conditions realistes. Le traitement mis 
en jeu consiste h extraire les donnees de variation de fre- 
quence due i I'effet Doppler, caracteristiques du degre de 
stenose, de 1* artefact cr£e par le mouvement de la. parol 
35 du vaisseau, le d£placement du volume de l'£chantillon, et 
les parasites de l f installation. 

Traditionnellement, le diagnostic de 1 f ar terippathie 
n4cessitait une arteriographie envahissante ( angiographic 
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de contraste). Plus recemment, la scanographie par ultrasons 
a conduit au diagnostic du mal sur la base d'affichages des 
spectres de transformation de Fourier analyses par un analyseur 
de spectres. Une telle lecture est difficile. Les spectres 
enregistres (tels qu'ils sont represents, par exemple, dans 
r. Knox, et al. , "Computer Based Classification of Carotid 
Arterial Disease : A Prospective Assessment," Stroke, Sept.- 
Oct. 1982, 13, 589-94 sont un elargissement continuel du 
signal a mesure que le degre de stenose augmente. Lire les 
spectres pour quantifier le degre de stenose exige une grande 
competence et une grande experience. Deceler un etat qui 
evolue lentement est presque impossible. 

Mais avec la methodologie de la presente invention, 
les donnees d" amplitude et de frequence de la transformation 
15 rapide de Fourier discrete sont traitees dans 1 • ordinateur 
a l'aide d ' algorithmes de reconnaissance des formes. Le pro- 
cede consiste a traiter les donnees des spectres enregistres 
en extrayant des morceaux utiles ou des caracteristiques 
predominantes des donnees. Avec le present systeme, le diag- 
20 nostic est obtenu objectivement et automatiquement a partir 
des donnees materielles., ce qui permet d'eliminer la neces- 
sity d'une grande competence ou d'une grande experience dans 
1 • interpretation . 

On technicien peut poser immediatement le diagnostic 
25 sans qu'un radiologue soit necessaire pour lire le spectre, 
ce qui entrainerait des retards de diagnostic. L' interpre- 
tation visuelle de caracteristiques telles que la frequence 
systolique de crete, la frequence diastolique, le degre d' elar- 
gissement spectral, ou la forme globale des formes d'ondes, 
30 rend virtuellement impossible, meme pour les radiologues 

les plus competents, une classification objective des di verses 
stenoses et la quantification des caracteristiques qui sont 
associees a chaque degre de stenose. La reconnaissance des 
formes des donnees d* amplitude et de frequence off re une 
35 methodologie statistique et analytique tres perf ectionnee 
permettant de dissequer les spectres Doppler a ultrasons 
pour en extraire automatiquement les donnees importantes 
pour le diagnostic. 
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Pour diagnostiquer automatiquement et objectivement 
le degre de stenose dans un vaisseau, dans l'artere carotide 
par exemple, un scanographe a ultrasons cree un signal Doppler 
a impulsions retrodiffuse en provenance du voisinage du flux 

5 sanguin central en plusieurs emplacements le long du vaisseau 
pres de l*entaille dicrotique. Les frequences de transfor- 
mation de Fourier de ce signal Doppler sont analyses dans 
un analyseur de spectres de transformation rapide de Fourier, 
en temps reel. Les spectres analyses en frequence sont en 

10 outre trait£s dans un ordinateur numerique specialise pour 

que la reconnaissance des formes entre une moyenne d' ensemble 
des spectres du patient et des spectres connus semblables 
et classes a l'avance dans une base de donnees puisse etre 
eff ectuee pour classer les donnees du patient dans une cate- 

15 gorie refietant le degr£ de stenose. Le systeme de classifi- 
cation est hierarchique, base sur des algorithmes de formation 
qui ponderent les mesures ou les "caractSristiques? des spec- 
tres. Pour que le traitement soit efficace, on comprime d'abord 
les donnees pour selectionner des parties ou des caracteris- 

20 tiques signif icatives du signal Doppler complet. La compression 
fait intervenir la creation de la moyenne d" ensemble, la 
selection d'un spectre ou d'un contour representatif d f ampli- 
tude en fonction du temps, a une frequence donnee representative 
du signal Doppler veritable, et 1 ' elargisseraent du spectre 

25 representatif avec des contours representatif s de frequence- 
amplitude autour du contour Doppler de base. 

La raise a la moyenne permet de r£duire l'effet des points 
de donnees ecartes, apparemment aberrants qui sinon pourraient 
poser un probl&me. Le cas ech^ant, on peut mettre a I'^chelle 

30 les largeurs des elargissements spectraux en divisant la 

largeur raesur£e par le cosinus de 1' angle Doppler (habituel- 
lement 60°). La mise a l'dchelle permet de determiner 1' inci- 
dence eventuelle de l'angle Doppler sur le diagnostic. On 
peut egalement explorer la sym^trie spectrale. 
35 La comparaison des caracteris tiques se fait avec le 

logiciel ARTHUR d' analyse de reconnaissance des formes, 
a 1' aides de methodes statistiques qui perraettent un diag- 
nostic empirique du degre de stenose. 
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Sur les dessins ci-joints : 

la figure 1 est un sch6ma-bloc simplifie de 1' instal- 
lation de diagnostic selon la pr^sente invention 

la figure 2 est une representation schema tique d*un 
signal Doppler type sur une courbe d •amplitude en fonction 
du temps ; et 

la figure 3 est une representation des donnSes comprim<5es 
dans une moyenne d' ensemble. Ces donn^es sont utilises dans 
le dispositif de reconnaissance des formes pour determiner 
le degre relatif de st6nose. 

On utilise un scanographe a ultrasons 10, de type Duplex 
pour creer la base de donn£es du patient. Le scanographe 10 
est un scanographe a effet Doppler combine en mode B et a 
impulsions, fonctionnant a 5 MHz, et fournit un signal Doppler 
en quadrature (bidirectionnel) . On peut supprimer la pratique 
courante consistant a filtrer le signal Doppler dans un filtre 
passe-haut pour Sliminer les parasites dus a la paroi. Les 
patients souffrant d'une art^riopathie ont souvent des compo- 
santes du signal Doppler proches de 0 Hz, et ces signaux 
seraient perdus avec un filtre passe-haut. Les signaux proches 
de 0 Hz sont prSdominants dans la p<5riode de fin de diastole 
ou au niveau de l'entaille dicrotique, oil une inversion de 
circulation peut parfois Stre observes. L ' electronique Doppler 
d'origine du scanographe commercial 10 avait son point 3 dB 
du filtre de paroi modifid de 1 kHz a 100 Hz pour que 1» infor- 
mation d'ecoulement ne soit pas perdue dans ces cas-la. On 
obtient ainsi une sensibility maximale aux signaux d'arterio- 
pathie en traitant une bande passante plus large. 

D' habitude, le scanographe 10 est un scanographe du 
type Mark V Duplex f ourni pax Advanced Technology Laboratories 
(ATL) , a Bellevue, Washington, Etats-Unis d'Am&rique (division 
des laboratoires Squibb). Ce scanographe utilise trois trans- 
ducteurs de 5 MHz a foyer fixe qui sont incorpores a une 
roue toumante, engendrant une image sector ielle bidimension- 
nelle des interfaces de tissus mous. L» image obtenue en mode B 
est utilisee pour identifier le vaisseau interessant, pour 
visualiser la calcification de la paroi, pour reconnaitre 
les variations anatomiques, pour placer le volume de I'echan- 
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tillon Doppler a impulsions au centre d'£coulement du vaisseau, 
et, ce qui est le plus important, pour maintenir un angle 
d' incidence Doppler type de 60° environ entre le faisceau 
Doppler et l'axe du vaisseau sanguin. L ' image obtenue en 

5 mode B est mise en memoire electroniqueroent, tandis que l'un 
des transducteurs est utilise dans le mode Doppler a impul- 
sions pour contrdler la vitesse du sang a l*int6rieur du 
volume d'^chantillon. L 1 orientation du faisceau Doppler et 
1 ■ emplacement du volume d'echantillon sont affich6s sur l 1 image 

10 obtenue en mode B. La tSte de balayage utilise habituellement 
des transducteurs de 6 mm de diametre mis au point a une 
profondeur de 25 mm pour la saisie des donn£es des art&res 
carotides primitive et interne. Le volume d'echantillon pour 
ce transducteur est de 3 mm 3 environ h cette distance focale, 

15 assurant une bonne limite de resolution au niveau de l»6coule- 
ment central oil les perturbations dues a la paroi ont peu 
d'effet. 

Bien que le dispositif soit utilisable pour la plupart 
des vaisseaux sanguins majeurs, les recherches s'attachent 

20 surtout actuellement au diagnostic du degr£ de st&iose dans 
les artferes carotides, et cette description sera limit^e 
a ces art&res. Les emplacements anatomiques a partir desquels 
des diagrammes de vitesse sont couramment enregistr£s pour 
!• analyse des art&res carotides sont I'artfere carotide primi- 

25 tive (1-CCA), l'artere carotide interne proche (p-ICA), et 
les emplacements situes le long de la bifurcation oil des 
signaux de haute frequence sont eventuellement d6tect6s. 
En chaque emplacement, 1' angle entre le faisceau Doppler 
et I'axe longitudinal du vaisseau est mesur6 directement 

30 a partir de 1' image obtenue en mode B. On obtient en mSme 
temps une information 12 de minutage d 1 electrocardiogramme 
(EGG) a l'aide d'un circuit qui se declenche a partir de 
l'onde R. En plus de la voix du technicien (un canal), des 
signaux de quadrature Doppler (deux canaux), ainsi que le 

35 signal d'onde R d'ECG module en frequence (un canal) sont 
enregistres automatiquement sur bande magn£tique a quatre 
canaux pour 1' analyse hors ligne. Environ 20 a 25 cycles 
cardiaques consecutifs sont enregistres sur bande a partir 
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de chaque emplacement. Quand une arteriopathie est detectee 
dans la partie £loignee de l'artere carotide primitive ou 
dans la partie proche de l'artere carotide interne {dans 
les deux premiers centimetres environ), c'est 1 1 emplacement 

5 de la perturbation maximale de l'ecoulement que l'on utilise 
pour 1' analyse. Ce signal est identifi6 par le technicien, 
soit d'apres une augmentation marquee des frequences Doppler, 
soit d'apres la "rudesse'* audible du signal* Les donnees 
enregistrees sur bande sont relues par 1 1 intermediate d'un 

10 analyseur de spectres en temps reel, et sa sortie numerique 
est en general ecrite .directement sur disque.. 

Dans les experiences de la presente invention, des donnees 
ont ete enregistrees a partir de l'ecoulement central de 
quatre emplacements pour chaque ensemble carotidien, donnant 

15 un total de huit emplacements par sujet : \es regions proche 
(a environ 3 cm du globe oculaire) et eloignee (a environ 
1 cm du globe oculaire) de l'artere carotide primitive, l'ar- 
tere carotide interne au niveau de l'ecoulement perturb^ 
audible maximal, et 1' emplacement de la carotide interne 

20 a 1 ou 2 cm de ce point. Le choix de deux emplacements de 
la carotide primitive etait base sur 1" experience qui indi- 
quait que des differences de comportement d'ecoulement pou- 
vaient etre constatees en differents points le long de cette 
artere. Sur la base d'essais complementaires, il a et£ deter- 

25 mine que les donnees provenant de l'artere carotide primitive 
proche et de l'artere carotide interne proche sont les seules 
donnees necessaires pour le dispositif prefer^ selon la pre- 
sente invention. 

La decomposition spectrale du signal Doppler directionnel 

30 est effectuee dans 1' analyseur de spectres 14 a 1'aide d'une 
methode de transformation de Fourier rapide discrete numerique 
(DFT). L'analyseur de spectres, tel que l'analyseur de spectres 
FFT en temps reel special qui est fourni par Honeywell, Inc. 
a Denver, Colorado, Etats-Unis d'Amerique, fournit 400 spectres 

35 par seconde avec une limite de resolution de frequence de 
100 Hz. Toutes les 2,5 msec, un nouveau spectre est cree, 
affichant 10 msec, de frequence Doppler. Les spectres adjacents 
comprennent ainsi 75 % de donnees anciennes et 25 % de donnees 
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nouvelles. L'intervalle d'analyse est de 10 kHz, 7 kHz" 6 tan t 
reserves pour l'ecoulement en sens direct et 3 kHz pour l'£cou- 
lement en sens inverse ou retour. L' existence de coraposants 
de frequence au-dela de 7 kHz se produit regulierement avec 

5 les stenoses severes, et peut etre deduite de la pente de 

l'enveloppe des donnees de la transformation de Fourier dis- 
crete (DPT) qui sont disponibles entre 0 et 7 kHz. D ' habitude , 
la pente de la DPT varie regulierement et se rapproche suf f i- 
samment de zero entre 0 et 7 kHz pour que le plus fort contenu 

10 de frequence des donnees Doppler soit Evident sur la base 
de 1* extrapolation des donnees de 0 a 7 kHz. 

Une interface de materiel et de logiciel avec un ordina- 
teur PDP 11/34 ou LSI 11/23 permet le transfert des donnees 
spectrales num^risees et de 1' information de minutage de 

15 l'ECG, directement de l'analyseur de spectre sur la memoire 

k disque 16. Un sous-programme de gestion de donnees du logiciel 
est congu pour permettre d'introduire dans le fichier de 
donnees les renseignements generaux concernant le patient 
(identification du patient , age, sympt8mes, bruits, operations 

2Q chirurgicales anterieures de la carotide, etc.), le cot£ 

et 1' emplacement du signal, I'angle Doppler et, eventuellement, 
d'autres informations. 

La decomposition spectrale du signal Doppler en quadra- 
ture se fait sous forme numerique. La precision de transfor- 

25 raation de ce signal est de 8 bits, et la configuration permet 
de delivrer 100 bits d 1 informations de frequence toutes les 
2,5 msec. La precision de I'analyseur de spectres est de 
+ 5 % en amplitude pour chaque bit de frequence. 

Le signal Doppler d ' ecouleraent sanguin affiche un echan- 

3q tillon statistique des mouvements des cellules sanguines 
individuelles a l'interieur du volume d'echantillon* Chaque 
pixel spectral, large de 100 Hz et long de 10 msec, repr6- 
sente une partie de cet echantillon statistique. La variation 
et la fluctuation d' intensity d'un pixel k 1' autre sont £tonna- • 

35 ment grandes. La variability statistique de chaque pixel 
se reflete dans la variability de deux spectres adjacents, 
meme s'ils partagent 75 % des donnees. 

Pour require la variance de 1 'estimation, on utilise 
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un systeme de mise a la moyenne spectrale. Ce systeme a 6t6 
coricu pour preserver la limite de resolution temporelle a 
I'interieur du cycle cardiaque tout en fournissant une meil- 
leure estimation de la repartition reelle des frequences 
Doppler. Le precede permet d'obtenir une famille de spectres 
mis a la moyenne, repr£sentant chacun une partie differente, 
de 10 msec, du cycle cardiaque. La mise a la moyenne renforce 
les caracteristiques de la forme d'ondes qui sont constantes 
ou periodiques, tout en r£duisant les autres caracteristiques, 
non periodiques ou variant avec le temps (variations de batte- 
ment a battement, parasites du dispositif non immobile, et 
autres artefacts). La mise a la moyenne permet d'£viter la 
prevention qui serait mise en jeu dans line selection subjec- 
tive d'un battement cardiaque "representatif p . 
15 On effectue habituellement la mise a la moyenne spectrale 

en introduisant les donnees de signal du disque dans l'ordi- 
nateur numerique specialise 18, au debit de donnees de 40 kilo- 
octets/sec . de l'analyseur de spectres, ou a tout autre debit 
de donnees convenable, en faisant la somme des signaux, et 
20 en divisant par le nombre total de signaux- On obtient 1* infor- 
mation de minutage de l'ECG a l'aide d'un circuit de materiel 
qui se declenche k partir de l'onde R, tout en rejetant les 
variations d' amplitude et de base et les ondes P et T accen- 
tuees. Les donnees de variation de frequence due a l'effet 
25 Doppler sont ainsi segment^es en cycles cardiaques correspon- 
dants bases sur le minutage de l'onde R. 

Etant donne que le retard electromecanique myocardique 
est relativement constant d'un battement a .1" autre, c'est 
l'onde R que l'on choisit comme reference de temps pour syn- 
chroniser la mise a la moyenne des donnees. Pour eiiminer 
les artefacts dus aux arythmies cardiaques, on rejette tout 
signal ayant un intervalle R-R qui ou bien diff&re de l'inter- 
valle precedent de 110 msec, environ en plus ou en moins, 
ou qui est plus long de plus de 33 % environ que 1' intervalle 
35 precedent. On pourrait utiliser d' autres facteurs comme r&gle 
de separation pour les arythmies. 

Les battements rejetes peuvent etre des battements 
arythmiques (irregularis cardiaques) ou anormaux, tels 
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que contractions ventriculaires prematurees. La r&gle des 
33 % comprend, cependant, tous les battements qui sont des 
variations normales du rythme cardiaque. Ici encore, on rSduit 
la prevention de 1» analyse, etant donne que l'on utilise 
une r&gle objective. 

On fait la somme des enregistrements spectraux obtenus 
a partir de plusieurs battements cardiaques par ailleurs 
consecutifs et acceptables, pour obtenir un spectre de moyenne 
d' ensemble pour le patient. On calcule la moyenne d' ensemble 
k l'aide de la formule : 

K 

E (f,t) = I B (f,t,n)/k 
n=l 

dans laquelle : E est l'amplitude DPT if, t, n ; 
f est la frequence ; 

t est le temps (mesure sur 280 intervalles 
de 2,5 msec, aprfes l'onde R de l'ECG) ; 

n est le numero d'ordre du spectre cardiaque ; 

et 

K est le nombre des spectres dans la moyenne 
d'ensemble (habituellement 15 i 20). 

Quinze a vingt spectres individuels fournissent une moyenne 
stable representative de l'^tat du patient. Si l'on fait 
la moyenne de moins de 15 spectres, il peut en rdsulter un 
diagnostic incorrect/ etant donne que la variance de la 
moyenne risque d'etre trop grande. Si l'on utilise plus de 
20 spectres, on constate une faible amelioration, et il peut 
en resulter une certaine reduction de la valeur du diagnostic. 
En effet, les spectres deviennent irr£guliers k cause de 
la fatigue du technicien qui deplace le scanographe, ou a 
cause de modifications du patient dues 4 la longueur de 
l'examen. Une moyenne de 20 spectres est pr6f Arable. 

On fait la moyenne de spectres acceptables consecutifs 
pour eviter d'introduire une distorsion due k la selection 
arbitraire des spectres. 

Une estimation du spectre de puissance reel du signal 
Doppler (obtenue en effectuant une transformation de Fourier 
discrete (DFT) ) est sujette a erreur statistique. On peut 
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etudier analytiquement le comportement de cette erreur en 
deduisant sa variance. Oppenheim et Schafer ont montre .que 
pour un proc<5de aleatoire gaussien : 

Var [f(f)l = P(f)[l + (sin 2 fN)/(N sin 2 f)| 
5 « P(f)|l + sine (2 fN) j 

sine (2 £) 

ou : P(f) est la densite spectrale de puissance estim£e 
(periodogramme) en fonction de la frequence f 
P(f ) est la density spectrale de puissance r<£elle (sous- 
10 jacente) a partir de laquelle les estimations observes 

|P(f)l sont engendrees 

N est le nombre de points de donnees dans la DFT 
sinc(x) = sin(x)/x. 
Le signal Doppler doit etre gaussien £tant donn£ qu'il est 
15 produit par un nombre assez grand de diffuseurs al^atoires 
(cellules sanguines). Comme : 
Lim Var lP(f ))1 = P(f ), 
. N » 

1' erreur de 1' estimation ne converge pas vers zero a mesure 

20 que le nombre d'^chanti lions, N, augmente (e'est-a-dire qu'il 
ne s'agit pas d'une estimation "coherente* ) . On peut cependant 
obtenir une amelioration en faisant la moyenne d • estimations 
successives. En particulier, Welch a montr6 que : 
Var (P(f)] = P(f)/K (si N est assez grand) 

25 oil : K est le nombre de moyennes spectrales. 

Comme la variance de 1' estimation, k chaque frequence, est 
egale au spectre de puissance r6el a cette frequence si l'on 
n'effectue aucune moyenne (K = 1), la DPT peut differer nota- 
blement du spectre sous-jacent r£el a moins que l'on utilise 

30 une moyenne. Une variance de la nature cons tat £e ici signifie 
que le spectre r£el est compris entre- 0 et deux f ois celui 
de 1* estimation qui est mesur^e. 

La reduction de variance fait la moyenne de spectres 
de 20 cycles cardiaques acceptables successifs, £ des instants 

35 identiques par rapport a chaque onde R. La taille du fichier 
de donnees resultant, mSme s'il ne represente que le 1/20 
des donnees brutes, est telle que la memoire a disque ne peut 
encore conserver de donndes que pour onze patients environ; 
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Une nouvelle compression de donnees est necessaire, et on 
va l'expliquer. 

Le progiciel ARTHUR, qui est fourni par 1 'University 
de Washington, est utilise comme base de sous-programmes 

5 pour effectuer de la reconnaissance de formes. Ce progiciel 
fournit le logiciel necessaire pour aborder le probleme gene- 
ral pose par la reconnaissance des formes dans des donnees 
stochastiques, lorsqu'on ne sait rien ou presque rien de 
la forme des densite*s de probability sous- jacentes. On met 

10 en oeuvre un systeme de gestion de base de donnees pour orga- 
niser les grandes quantity s de donnees qui sont raises en 

L'approche reconnaissance des formes est basee sur 1' uti- 
lisation de donnees representant des etats connus de la nature 

15 (maladie par exemple) pour "entrainer" un systeme de classi- 
fication. On utilise ensuite cette "base de donnees de forma- 
tion" pour mettre au point une regie de decision, que l'on 
utilise pour classer des donnees inconnues. 

La base de donnees de formation, qui est fournie par 

20 1' University de Washington et la Washington Research 

Foundation, se compose de donnees provenant d'arteres malades 
et normales. Des patients considers comme souffrant d'arte*- 
riopathie ont yte" evalues independamment par arteriographie 
de contraste ( angiographic ) , qui a fourni une estimation 

25 du pourcentage de reduction de diametre & partir de la moyenne 
de mesures independantes prises par deux radiologues mais 
lues par un seul d'entre eux. L 1 estimation utilisee etait s 
% de stenose = 100 X (1. - diametre du vaisseau non bouchyj 

diametre total du vaisseau / 

30 L'age des patients malades allait de 38 & 82 ans, avec une 
mediane de 61 (moyenne = 61,7). On a obtenu des donnees nor- 
males a partir de sujets qui dtaient presumes non malades 
en raison de 1* absence de symptSmes et de bruits J Leur 3ge 
s'^chelonnait de 21 a 43 ans, avec un Sge median de 27 

35 (moyenne = 28,5). La base de donnees de formation comprend 
des donnees comprimees de spectres de moyenne d' ensemble 
de chaque etat connu et les donnees du signal Doppler k ultra- 
sons complets sont disponibles pour une validation croisee 
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ou une etude £largie, si n^cessaire. 

La base de donnees de formation ne sera probablement 
jamais reassemble etant donnd que l'arteriographie de con- 
traste n'est pas regulifcrement prescrite en raison des risques 
gu'elle presente. Avec la supr&natie des techniques de diag- 
nostic par ultrasons et autres techniques hon envahissantes, 
les medecins se refusent & autoriser une art£riographie de 
contrast e qui impose aux patients des risques semblables 
£ une operation chirurgicale. 

Les etats pathologiques qui sont identif iables k l r int<5- 
rieur des donnees Doppler sont consid6r£s comme cat£goriels, 
plutot que continus. S'ils sont categoriels, l f approche 
de diagnostic : (1) reconnait que des groupements naturels 
peuvent exister dans les donnees (par exemple, normal ou 
malade) r (2) permet une estimation de la probability d'erreiir 
de classification, et (3) r£duit le problfeme d'erreur dans 
les lectures individuelles des arteriogrammes en cherchant 
un groupage naturel des spectres correspondant a un 6tat 
particulier des donnees, plut3t qu ' en . cherchant une varia- 
tion pr<*visible d'une facon regulifere k mesure que la st&iose 
augmente. Des erreurs peuvent r^sulter, dans les spectres 
individuels, d' incertitudes inh^rentes sur le diam&tre du 
vaisseau non bouch£, meme avec la mesure par artSriographie 
de contraste. Les estimations de la stenose sont bashes sur 
des images planes d'une lesion tridimensionnelle irrdguli&re 
£ I'interieur du vaisseau. Une etude r6cente indique que 
des taux de disaccord entre radiologues de 17,5 % peuvent 
se produire quand on utilise 1 • art^riogramme pour identifier 
un etat pathologique comme une stenose h 0 %, 1-9 %, 10-49 %, 
50-99 %, ou 100 %, ou de 12,5 % quand on identif ie un etat 
pathologique comme une stenose h 0 %, 1-20 % ou 21-99 %. 
Ces variations sont r£duites avec une m£thodologie cat<Sgo- 
rielle de reconnaissance des formes > qui cherche tout simple- 
ment un groupage etant donnS que les tendances majeures sont 
reconnaissables meme si de fins details peuvent dchapper 
& 1 1 observation . 

Les donnees du patient sont classees selon une serie 
hierarchique de decisions binaires separSes, selon l'un des 
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sch&nas suivants : 

A 

Pour un vaisseau qui n'est pas complfeteraent bouche. 
Decider s'il est normal ou malade ; 
5 S'il est malade , decider si la st^nose est supdrieure 

ou inferieure a 50 % ; 

Si la stenose est inferieure & 50 %, decider si elle 
est supdrieure ou inferieure k 20 %. 

B 

W Pour un vaisseau qui n'est pas completement bouche, 

Decider s'il est normal ou malade (partiellement bouche) ; 
S'il est malade, decider si I'occlusion est superieure 
ou inferieure a 50 % ; 

Si l f occlusion est inferieure a 50 %, decider si l'occlu- 

15 sion est superieure ou inferieure a 35 % ; 

Si !■ occlusion est inferieure & 35 %, decider si 1» occlu- 
sion est superieure ou inferieure a 15 %. 

Le cas echeant, on peut realiser une division en cate- 
gories plus fines. 

20 La decision initiale doit etre : le vaisseau est-il 

ouvert ou completement bouche (occlus) ?, etant donne qu'une 
erreur de diagnostic forc£e se produira si le vaisseau est 
bouchS et ne laisse apparaitre aucun ecoulement. La hierarchie 
de la decision ulterieure presume un ecoulement dans le vais- 

25 seau et impose les decisions ulterieures. On procede a la 
determination initiale (ouvert ou bouche) en observant un 
signal d 1 ecoulement dans l'art&re carotide interne ainsi 
que dans l'art&re carotide primitive proche. Si un ecoulement 
est detecte dans la carotide interne, on peut utiliser pour 

30 le diagnostic le dispositif informatique de reconnaissance 
des formes. Si aucun ecoulement n'est detecte, c'est que 
le vaisseau est complfetement bouche, et le technicien pose 
le diagnostic sans utiliser la reconnaissance des formes 
spectrales. 

35 Quand on comprime les donnees mises a la moyenne, on 

peut extraire les informations essentielles necessaires pour 
une analyse reussie de reconnaissance des formes. Ces "infor- 
mations essentielles" se composent d 1 estimateurs d 1 emplacement 
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et d'echelle pour les spectres Doppler, qui sont des indica- 
teurs d'une frequence a l'interieur du signal et de l'etale- 
ment du spectre associ6 a cette frequence caractexistique* * 
Pour en tamer la compression de donnees, il faut trouver 
dans les parasites regus avec le signal rdtrodif f use, un 
signal Doppler representatif de l'ecoulement sanguin (et 
indicatif du degr£ de stenose). Un algorithme permettant 
de remplir cette fonction a 6t6 mis au point et a l'epreuve 
sur des donnees cliniques. Cet algorithme permet de trouver 
la frequence de mode (amplitude maximale) du signal Doppler, 
apres utilisation d'une statistique d'ordre croissant ou 
decroissant pour trouver la region du spectre conteriant le 
signal. La statistique d'ordre croissant ou decroissant 
resiste a des parasites a bande passante dtroite de grande 
amplitude (artefact de paroi par exeraple). La frequence d 1 am- 
plitude maximale pour la moyenne spectrale determine la fre- 
quence de mode de signal permettant d'extraire des 280 inter- 
valles de temps de la moyenne d' ensemble le signal de contour 
caracteristique d 1 amplitude en fonction du temps. Autrement 
dit, pour chaque interval lede temps dans le cycle cardiaque 
enregistre" dans les donnees de variation de frequence Doppler, 
on utilise 1' amplitude de la frequence d'amplitude maximale 
pour definir le contour de signal de mode en trouvant 1* ampli- 
tude du signal a cette frequence dans chacun des 280 inter- 
valles de temps (pixels) d'un bout a 1" autre du spectre; 

On determine la frequence d f amplitude maximale en tragant 
une courbe de 1' amplitude en fonction de la frequence a partir 
des donnees de variation de frequence Doppler-. En tant que 
telles, 1' amplitude maximale et sa frequence correspondante. 
sont faciles a extraire. 

II est commode d'estimer l'echelle ou la largeur du 
signal d'ecoulement de variation de frequence Doppler pour 
le contour caracteristique en ajustant des contours a 3 db 
(les points de demi-puissance du signal) et a 9 db du contour 
de signal de mode. On ameliore les estimations du contour 
de signal de mode et du contour d'etalement en utilisant 
dans le calcul un filtre de mecliane courante. Cette approche 
s'est revelee insensible a la plupart des cas d r artefacts 
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de paroi, de bruits de fond, et d' interferences de canal 

en quadrature. Bien entendu, on peut utiliser d'autres tech^ 

niques de mise a l'^chelle ou d'elargissement. 

Etant donne que 1" information d'ecoulement en fin de 

5 diastole est susceptible d'etre interessante dans l'analyse, 
on a besoin d'une m£thode d'obtention de cette information 
a partir du spectre de moyenne d f ensemble. La methode doit 
etre insensible aux intervalles R-R, variant largement, qui 
sont observes entre des patients differents. Mais, en incluant 

10 100 msec, de donnees priced ant immediatement chaque spectre 
d'onde R acceptable dans la moyenne d* ensemble, on peut aise- 
ment inclure 1 • information de fin de diastole, dans les con- 
tours comprimes representatif s de la bande passante de fre- 
quence et d' amplitude des donnees de variation de frequence 

15 Doppler completes. 

Les donnees comprimees (fig. 3) comprennent cinq lignes 
de contour representant 1* amplitude Doppler en fonction du 
temps pour chaque frequence. La ligne de contour centrale. 100 
(fig. 3) represente la frequence mediane ou frequence de 

20 signal de mode qui represente le mieux I'ecoulement sanguin 
au centre du vaisseau. La paire de lignes ou de "contours 
d'£chelle M 200 ou 300 entourant immediatement le contour 
de signal de mode 100 signale les points h 2 db, et la 
paire ext^rieure extreme 400 et 500 represente les points 

25 ^ 9 db par rapport au mode. Ces contours -compriment les don- 
nees a 1400 octets par emplacement (f 4. contours + 1 mode) x 
(280 intervalles de temps 1 = 1400 octets). On obtient une 
nouvelle reduction de donnees de vingt k un. Les donnees 
comprimees pour plusieurs centaines de patients peuvent 

30 s' adapter sur un seul disque, ce qui facilite l f analyse des 
donnees. 

Le systfeme de reconnaissance des formes permet de compa- 
rer des caracteristiques de trois types generaux : (1) celles 
qui font intervenir des relations avec le temps de la forme 
35 d'onde de vitesse ; (2) celles qui font intervenir la largeur 
du spectre en des points choisis a l'avance le long du cycle 
cardiaque ; et (3) celles qui font intervenir une decomposi- 
tion de frequence de la forme d'onde de frequence moyenne 
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estimee. La forme d'onde de frequence moyenne ou la ligne 
de contour de mode est analysee dans une transformation de 
Fourier discrete a fenStre pour extraire les composantes 
de frequence Doppler principales entre 1 et 20 Hz. Quatre- 
5 vingt-quatorze caracteristiques ont ete initialement iden- 
tifies comroe des caracteristiques potentielles pour compa- 
raison, et sont decrites dans le tableau I. 

TABLEAU 1 

Caracteristiques potentiellement imoortantes de 
" l0 reconnaissance des formes 

On a mis a l'echelle ou transforms certaines de ces 
caracteristiques, a l'aide de fonctions logarithmiques ou 
de puissance, pour obtenir les caracteristiques utilisees 
effectivement dans les equations discriminantes. 
15 I « Artere carotide primitive proche 

(1) Aire sous la forme d'onde d'ecoulement moyenne 

(2) Cos (e) 

(3) Sin (0) 

(4) Aire sous le pic systolique/cos (9) ± intervalles 
de 15 msec. 

A - Caracter istiques de frequence de mode 

(5) Pin de diastole (50 msec, avant l'onde R de l'ECG) 

(6) Pic systolique 

(7) Premier point de pente zero apres le pic systolique 

(8) Encoche dicrotique 

(9) Diminution de frequence du pic systolique au premier 
point de pente zero apres le pic systolique/cos (e) 

B « Caracter istiques de debit diastolique 

(10) 125 msec, apr^s 1' encoche dicrotique/cos (e) 

(11) 50 msec, apres l'onde R/cos (6) 

(12) Frequence d' amplitude maximale dans la transforma- 
tion de Fourier discrete 4 fenetre (WDPT) entre 
harmoniques 4 et 9 

(13) Frequence d' amplitude maximale dans la WDFT entre 
les harmoniques 5 et 9 

(14) Frequence d' amplitude maximale dans la WDFT entre 
les harmoniques 6 et 9 

(15) Rapport de Pourcelot 
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(16) Frequence systolique de mode/frequence de d£but 
de diastole de mode 
c - Caract6r istiques reflgtant la presente ou 1' absence 
d'une "fenetre systolique" 
5 (18) Frequence maximale inf£rieure a 3 db entre t 

(systole - 30 msec.) et t = (systole + 50 msec.) 
(19) Frequence .maximale inferieure a 9 db entre t = 

(systole - 30 msec.) et t = (systole + 50 msec). 
D. Caractgristiques d' information de phase (annoncfees) 
10 (20) Phase du second harmonique 

(21) Phase du troisifeme harmonique 

(22) Phase du quatrieme harmonique 

E * Caractfer istiques basees sur un elargissement spectral 
(Pour les contours de ± 3 db et ±9 db f respectivement ) 
15 (23-26) Fin de diastole 

(27-30) Systole 

(31-34) Premier point de pente zero apres le pic systo- 
lique 

(35-35). 100 msec, aprfes la systole 
20 (39-42) Systole/cos (e) 

(43-46) 75 msec, aprfes la systole/cos (0) 
(47-50) 100 msec, aprfes la systole/cos (e). 
II. Caract6ristiques de la carotide interne proche 

(I) Aire sous la forme d'onde d'ecoulement moyenne 
25 (2) Cos G 

(3) sin 6 

(4) Aire sous le pic systolique/cos (e ) ± intervalles de 
15 msec. 

A. Caracteristiques de valeur de mode 

30 (5) Fin de diastole (50 msec, avant l'onde R de l'ECG) 

(6) Pic systolique 

B. Caracteristiques maximales 

(7) Frequence 3 db de crSte 

(8) Frequence 9 db de crSte 

35 (9) Degre de deceleration suivant la systole/cos (e) 

(10) Energie de 1' amplitude maximale entre les harmo- 
niques 5 et 9 

(II) Frequence d' amplitude minimale dans la WDFT de 2 
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au maximum precedent 
(12) Rapport de Pourcelot = (systole - ler point de 

pente zero)/raoyenne) , en utilisant la forme d'onde 
de mode 

5 (13) Rapport = (systole - ler point de pente zero ) /systole ) , 

en utilisant la forme d'onde de contour de 9 dB 
superieure 

(14) Rapport = (systole, en utilisant la forme d'onde 
moyenne estimee)/fin de diastole, en utilisant 

10 la forme d'onde de 9 dB superieure 

C Caracteristiques refletant la presence ou 1 'absence 
d'une "fenetre systolique" 

(15) Frequence maximale 3 db infSrieure entre t = 
(systole - 30 msec.) et t = (systole + 50 msec.) 

15 (16) Frequence maximale 9 db infSrieure entre t « 

(systole - 30 msec.) et t = (systole + 50 msec:) 
D. Caracteristiques bashes sur un glargissement spectral 
(Pour les contours de ±3dB et ± 9dB f respectivement ) . 
(17-20) Fin de diastole (50 msec, avant l'onde R) 
20 (21-24) Pic systolique 

(25-28) 50 msec, apres le pic systolique 
(29-32) 100 msec, aprfes le pic systolique. 
(33-36) Pic systolique/cos (e) 

(41-41) 100 msec, aprfes le pic systolique/cos ( e) 
25 ou e est 1' angle forme par le faisceau Doppler avec l'axe 
du vaisseau. 

Pour les caracteristiques faisant intervenir des rela- 
tions avec le temps, on a encore comprint le spectre de 
moyenne de 1' ensemble de 5 contours en un seul contour moyen 
30 pondere, les poids etant affectes en fonction des amplitudes 
relatives du contour en question, seion la formule : 
f(t) - (m(t) + 2 [Ui (t) + Li (t)J + 

b(U 2 (t) + L 2 (t)J}/ (1 + 2 [a + bj) 
dans laquelle : 
35 m(t) est le mode spectral au temps t ; 

Ui(t) est la frequence 3 dB superieure au temps t ; 
I*i (t) est la frequence 3 db inferieure au temps t ; 
□2^) est la frequence 9 dB superieure au temps t ; 
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I*2(t) est la frequence 9 dB inferieure au temps t ; 
a = 10-3/20 . et 
b = 10-9/20. 

On a situ£ le pic systolique et la premiere inflexion du 
5 pic systolique a l'aide de la premiere derivee lissee mediane 

du contour 100 de frequence (ou de mode) d 1 amplitude maximale. 

On a lisse ce contour selon la methode de Mosteller et Tukey. 
Pour les caracteristiques faisant intervenir la largeur 

du spectre, on a mesure la largeur du contour de mode au 
10 contour d'echelle voulu et on en a fait la moyenne pour une 

fenetre de ± 12 , 5 msec. (5 intervalles) centree sur le temps 

specific. 

Comme le dispositif de reconnaissance des formes pondere 
les caracteristiques de facon lineaire, il est important 

15 de transformer de facon non lineaire les "caracteristiques 
brutes" (celles qui sont choisies dans le tableau I) pour 
tenter de porter au maximum leur efficacite dans les equations • 
discriminantes lineaires. Par consequent, on a calcule des 
fonctions exponentielles et de puissance des caracteristiques 

20 de base, et on a trace les courbes correspondantes pour veri- 
fier leur linearite. Pour la commodity et la siraplicite de 
la presente invention, on a utilise une regression lineaire 
pas k pas pour selectionner les caracteristiques, et cette 
methode s'est r^velee adequate dans la mesure ofc les carac- 

25 t^ristiques selectionnees contiennent des informations uniques 
et non des informations redondantes. 

On a remis a I'echelle les caracteristiques pour obtenir ■ 
une variance nulle de la moyenne et des unites. Par analyse 
discriminante pas & pas des 94 caracteristiques, on a op6r6 

30 une classification des vaisseaux connus dans la base de don- 
nees (ensemble de formation) en utilisant la meilleure combi- 
naison de 5 caracteristiques selectionnees parmi les 94 pour 
chaque decision hierarchique. Ensuite, on n'utilise que ces 
5 caracteristiques pour diagnostiquer le degre de stenose 

35 £ parti r de la moyenne d 1 ensemble d'un patient. Les 5 carac- 
teristiques utilises pour chaque decision hierarchique (a 
savoir normal ou malade ; stenose superieure ou inferieure 
a 50 %, etc.) sont selectionnees separement par I'analyse 
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discriminante parmi 1' ensemble des 94 caracteristiques d'ori- 
gine, par rapport aux donn^es des etats art^riens connus. 
Bien entendu, on pourrait utiliser d'autres caracteristiques 
ou des caracteristiques supplementaires, s^lectionnees de 

5 fa9on differente. On pr^fere utiliser 5 caracteristiques 
car elles semblent fournir un haut degre de precision. On 
obtient un equilibre assez acceptable pour que chaque carac- 
teristique soit significative, sans que soient introduites 
des caracteristiques relativeroent peu signif icatives ou redon- 

10 dantes. Bien que la methode puisse etre orients dans un 
sens ou dans 1" autre par la selection de caracteristiques 
sur la base de cette analyse discriminante, 1' Evidence empi- 
rique montre que ce format repr£sente un outil extremement 
precis, automatique et objectif pour diagnostiquer correcte- 

15 ment la st^nose des artferes carotides. 

Les caracteristiques pr^ferees seiectionnees pour chaque 
decision sont : 

TABLEAU II 

1. Normal ou malade 
20 . In (de la valeur de l'ecoulement moyen de debut de 

diastole dans l'artire carotide primitive/cos e ) 
. - [variation post-systolique de frequence/cos el 4 au 
niveau de l'artfere carotide primitive 

. In (mode de la forme d'onde DPT) au niveau de l'artfere 
25 carotide interne 

. (point 9 dB inferieur de la forme d'onde} 4 au niveau 
de l"art&re carotide interne 

. In (ecoulement moyen de fin de diastole/cos 0 ) au niveau 
de l'art&re carotide primitive 
30 2. Stenose superieure ou inferieure a 50 % 

. f frequence 9 dB de crete] 2 au niveau de I'artere caro- 
tide interne 

. [frequence 9 dB inferieure maximale entre ±50 msec, 
du pic systoliquel 4 au niveau de l'artfere carotide 
35 interne 

. -[variation post-systolique de frequence/cos e J au 
niveau de I'artere carotide interne 

. In [premier minimum de la forme d'onde DPT I au niveau 
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de l'art&re carotide interne 

. Ilargeur 9 dB sup^rieure 50 msec, apr&s la systole l 4 
au niveau de l'artere carotide interne 
Stenose sup^rieure ou inf£rieure & 35 % 

• frequence 9 dB de crgte au niveau de I'art&re carotide 
interne 

. In I frequence 125 msec, avant I'encoche dicrotique 
de I'artere carotide primitive! 

. In (variation de la frequence systolique au niveau 
de l'art&re carotide interne/cos 01 
. In I forme d'onde DFT 9 dB inf^rieurej au niveau de 
l'art&re carotide interne 
Stenose superieure ou inf^rieure a 15. % 

. In (aire situee sous le pic pour la periode ±37,5 msec, 
autour de la systole I au niveau de l'artfere carotide 
primitive 

. In [frequence 9 dB de crSteJ au niveau de l'artfere 
carotide interne 

. frequence 9 dB de crSte ail niveau de l'art&re carotide 
interne 

. largeur 9 dB superieure 125 msec, aprfes la diastole 
au niveau de l'artfere carotide primitive 
. In (largeur 3 dB inf^rieure 125 msec, aprfes la diastole] 
au niveau de l'artfere carotide interne 
Stenose superieure ou inf£rieure a 20 % 

. In [aire situ£e sous le pic pour la periode ±35,5 msec, 
autour de la systole 1 au niveau de l'art&re carotide 
primitive 

. In (frequence 9 dB de cr§te| au. niveau de I'art&re 
carotide primitive 

. frequence 9 dB de crete au niveau de l'artfere carotide 
interne 

. largeur 9 dB inffirieure 100 msec, apres la systole 
au niveau de l'artfere carotide primitive 
. In (mode de la forme d'onde DFTJ au niveau de I'art&re 
carotide primitive. 

Cette determination de 20 % de stenose -est apparemment 
ez precise si l'on utilise un essai redondant en ce qui 
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concerne la frequence 9 dB de crete. 

S'il est souhaitable de determiner si la stenose vascu- 
laire est super ieure ou inferieure a 80 %, les caract£ris- - 
tiques preferees sont : 

6. Stenose super ieure ou inferieure a 80 % 

. In llargeur 9 dB superieure, a la systole] au niveau 
de l'artere carotide interne 

» In [frequence' 3 dB maximale entre ± 50 msec, du pic 
systolique] au niveau de l'artere carotide interne 
. In [premier moment du pic systolique/f requence 9 db 
superieure 125 msec, apres la diastole] au niveau de 
l'artere carotide primitive * 

. In (mode de la forme d'onde DFT] au niveau de l'artere 
carotide interne. 
15 On utilise une methode de resolution hyperplane de 

la reconnaissance des formes bas£e sur 1" analyse de regression 
lineaire dans le logiciel ARTHUR. La methode pr£f£r£e selon 
Pietrantonio et Jurs (generaleraent connue sous le nom de 
Pjurs) consiste a minimiser de facon iterative la fpnction : 

|Y - tanh(Wo + 2 WiXj)] 
dans laquelle 

Y « +1 si Ton sait que les donnees se trouvent dans 

la categorie 2 
=* -1, si l'on sait que les donnees se trouvent dans 
25 la categorie 1 

Y = categorie predite 

K = nombre de caracteristiques 
= poids de la caracteristique 
a valeur mesuree de la caracteristique i 
w Q = const ante. 

Pour s 'assurer un plus haut degre de precision/ on uti- 
lise une validation crois^e. La validation crois^e consiste 
a diviser la base de donnees de formation en deux sous- 
ensembles pour chaque categorie de stenose experimented, 
a "entrainer" ou former le classif icateur sur la base d'un 
sous-ensemble, puis a r£essayer la regie de decision du clas- 
sificateur en utilisant le classif icateur sur le second sous- 
ensemble de donnees. On a choisi un nombre choisi au hasard 
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des donnees spectrales comprim^es, comprenant 50 % environ 
de la base de donnees, a chaque iteration du processus de 
validation croisee, dans un syst&me de validation croise 
dit "dani-Schantillan plus ccrplement". On a restreint les sous- * 

5 ensembles pour eviter de surestimer la veritable precision 

de la r&gle de classification, comme cela pourrait se produire 
si le sous-ensemble choisi reprdsentait une partie plus iropor- 
tante de la base de donnees. 

La precision estimee du classif icateur obtenu avec ce 

10 precede de validation croisee est le pourcentage moyen correct 
dans les ensembles d'essai. A chaque decision dans la hi£rar- 
chie du format A, le pourcentage correct etait de 97 % pous 
la premifere decision (normal ou malade), de 95 % pour la 
deuxifeme decision (stenose superieure 6u inferieure & 50 %), 

15 et de 84/3 % pour la troisieme decision (stenose superieure 
ou inferieure & 20 %). On a utilise la somme des 5 caracte- 
ristiques ponder^es pour determiner le repfere discriminant 
de chaque decision. 

Le syst&me de classification actuellement utilise passe 

20 par quatre decisions binaires separees dans un crdre hi6rar- 
chique> a savoir : 

(1) Decider si le vaisseau est normal ou malade ; 

(2) s'il est malade, decider si la stenose est supe- 
rieure ou inferieure k 50 % ; 

25 (3) si la stenose est inferieure a 50 %/ decider si 

elle est superieure ou inferieure h 35 % ? et 

(4) si la stenose est inferieure a 35 I, decider si 
elle est superieure ou inferieure a 15 %. 

On etablit de la sorte six categories, k savoir : 

30 (1) Vaisseau non malade, (2) stenose de 0 a 15 %, (3) stenose 
de 16 a 35 I, (4) stenose de 36 a 50 %, (5) stenose de plus 
de 50 %, et (6) vaisseau corapl&tement bouche (stenose de 
100 %). Ces categories sont arbitraires pour les divisions 
intermediaires, et ont ete choisies sur la base de l'etendue 

35 des donnees obtenues a partir d'arteriographes, et du desir 
d'obtenir uhe finesse relativement grande d' evaluation de 
la progression de la maladie. On pourrait etablir d'autres points 
de rupture, tel que le format 0-20, 20-50 et plus de 50 %, 



25 



2591884 



dej& expose. La demarcation entre une st^nose sup^rieure 
ou inferieure a 15 % donne un avertisseraent prdcoce de la 
progression dangereuse de la maladie. Le radiologue a pose 
le premier diagnostic evalu6 a 15. % d' occlusion ou moins, 
aussi ce point est-il la premiere etape de la progression 
du mal qui peut etre obtenue avec cette base de donn£es. 

Mais 1'essai ultime des performances du classif icateur 
consiste a ^valuer les donn^es de patients souffrant d'une 
stenose connue, mais dont les donn^es n'ont pas 6t& utilises 
dans la conception du classif icateur . Les patients d'essai 
analyses pour essayer le classif icateur se composaient de 
neuf volontaires, qui 6taient presumes normaux / et de 95 pa- 
tients atteints d'une arteriopathie extracranienne evalu6e 
par angiographies Dans 1' ensemble du groupe, il y avait 
170 cotes convenant a l'£tude. Les 38 autres cotes n'ont 
pas ete inclus dans 1 'analyse pour les raisons suivantes : 
(1) 16 cotes presentaient une occlusion de 1'artfere carotide 
interne ; (2) 12 c6tes avaient subi une endart^riectomie avant 
1» etude aux ultrasons ; et (3) 10 cdt6s ont 6t& perdus a 
cause d'une defaillance mecanique du dispositif d'enregis- 
trement des donnees. 

On a proc£de k un essai en double aveugle dans lequel 
le radiologue, qui lisait les arteriogrammes de contraste 
biplans, n'etait pas au courant du diagnostic fourni par 
1' analyse inf ormatique , et le technicien Doppler n'etait 
pas au courant du diagnostic du radiologue. 

On a calcule la statistique Kappa telle qu'elle est 
utilisee par Cohen, Fleiss et Langlois, pour mesurer le degr£ 
d' accord entre deux specialistes evaluant ind^pendamraent 
la base de donnees. Compte tenu du hasard dans le diagnostic 
ou la decision, la formule devient ; 

Kappa = (p 0 - P e )/(1 - P e > 
oh P 0 est la probability d' accord observee et p e est 1' accord 
qui interviendrait sur la base du seul hasard. La quantity 
1 - P e me sure le degr6 d' accord qui peut etre obtenu par 
rapport et au-deli de celui qui serait pr6dit du simple 
fait du hasard. Le degre d' accord ef f ectivement atteint en 
plus du hasard est egal a P Q - P e . Si 1' accord observe est 
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supSrieur ou egal a l 1 accord resultant du hasard, Kappa sera 
superieur ou egal a zero, une valeur maximale de + 1 corres- 
pondant a un accord parfait. De meme, si 1" accord observe 
est inferieur a I'accord resultat du hasard, la valeur de 

5 . Kappa sera negative. Aux fins de ce calcul, on a suppose 
que la taille de 1 'echantillon 6tait grande. D'apr&s les 
r£sultats erapiriques, cette hypothese semble acceptable. 

L 1 accord entre les resultats de I'ordinateur et l'angio- 
graphie etait de 96,5 % (164/170) pour la decision normal 

10 ou malade, de 93 % (158/170) pour la decision stenose sup£- 
rieure ou inf6rieure k 50 %, et de 93 % (158/170) pour la 
decision st&nose superieure ou inf£rieure h 20 %. Des 29 
cot^s norraaux, 27 (93 %) ont ete correctement classes. Les 
deux autres cotes ont ete classes dans la categorie 1-20 %, 

15 et provenaient tous deux des sujets presumes normaux (ag£s 
de 17 et 34 ans), qui en fait n'avaient jamais subi d 1 angio- 
graphic de contraste. Par consequent, ils pouvaient effecti- 
vement souffrir d*une stenose mineure. 

Des 54 cotes classes dans la categorie 10 a 20 % de 

20 stenose, 44 (81,5 %) ont 6t6 correctement classes. Des 10 
cot£s restants, 4 ont ete classes comme normaux et 6 dans 
la categorie 21-50 %. Des 37 cotes classes dans la categorie 
21-50 %, 29 (78 %) ont ete correctement classes. Des 8 cotes 
restants, 5 ont £te mal classes dans la categorie 1 a 20 %, 

25 et 3 dans la categorie 51 & 99 % (pour tous les trois, on 

avait diagnostique par angiographic une reduction de diam^tre 
de 50 %). Des 50 cotes classes dans la categorie 51 a 99 %, 
41 (82 %) ont £te correctement classes. Des 9 cdt£s restants, 
1 a £te mal classe dans la categorie 1 a 20 %, et 8 dans 

30 la categorie 21 a 50 %. 

La precision globale de la classification informatique 
par catdgorie £tait de 83 % (141/170). La sensibilite de 
la m£thode.ou la faculte de reconnaitre la presence de la 
maladie 6tait de 97 % (137/141) et la specif icite de 93 % 

35 (27/29). Le degre d f accord entre la reconnaissance des formes 
et 1' angiographic, mesur6 par la statistique Kappa, £tait : - 
K = 0,769 ± 0,039 (erreur type). 

Seize cotes, atteints d* occlusion d'apr&s 1' angiographic. 
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ont 6te omis. Le diagnostic d'une occlusion a et6 reconnu 
par 1' absence de signal dans I'artere carotide interne, excluant 
ainsi une analyse par ordinateur. En aucun cas, un signal 
n'a ete soumis a 1' analyse par ordinateur quand 1'art^re 
5 avait ete dite occluse d'aprfes 1 ■ angiographic , etant donn6 

que le classif icateur ne pouvait que trouver un certain 6cou- 
lement dans le vaisseau. Mais aucun vaisseau n'a et6 omis 
dans l 1 analyse par ordinateur du simple fait qu'une occlusion 
avait et£ diagnostiqu^e au scanographe Duplex mais pas par . 
10 angiographic. 

On peut eventuellement a jouter a la hierarchie de. deci- 
sion la decision vaisseau complfetement bouche ou vaisseau 
ouvert, en effectuant la recherche par ordinateur pour un 
signal de variation de frequence Doppler dans I'art&re caro- 
15 tide interne. De la sorte, toutes les decisions seraient 

prises par 1 1 ordinateur , et aucun jugement ne serait laiss£ 
au technicien. 

Comme la variation de frequence due k I'effet Doppler 
est directement proportionnelle au cosinus de 1' angle Doppler 

20 par rapport a l'^coulement sanguin, on ne peut proc6der a 
un calcul de la vitesse d'ecoulement que lorsque cet angle 
Doppler est connu. La calcification de la. paroi et des varia- 
tions anatomiques (dont on sait qu'elles modifient toutes 
deux le signal d'£coulement) sont faciles a ^valuer par l'ima- 

25 gerie en mode B, 

Un appareil Doppler ^ impulsions, par d^clenchement 
periodique entre bornes, permet d'interrbger l'ecoulement 
a des profondeurs specif iques le long de I'axe du faisceau. 
L 1 experience clinique anterieure a montre qu'un scanographe 

30 ayant un point focal situe a environ 2 cm au-dessous de la 
surface de la peau convient pour £tudier les art&res caro- 
tides. D ' autre part f un petit volume d * echantillon par rapport 
au diametre du vaisseau permet l'^chantillonnage de l'ecoule- 
ment k partir d'un emplacement central, ce qui permet d'eviter 

35 les gradients de vitesse pr&s de la paroi du vaisseau. Le 
scanographe avait une distance focale de 25 mm, et on utili- 
sait un volume d'echantillon de 3 nmv* a son foyer pour obte- 
nir la limite de resolution voulue. L* image du scanographe 
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permettait d'etablir des reperes anatomiques utilises pour 
la normalisation des emplacements d'echantillonnage. Ces 
reperes etaient essentiels pour comparer les etats de patients 
diff brents ou pour comparer des etudes successives sur le 

5 meme patient. 

Etant donne que dans le procede de classification pre- 
£&r6 il n'a pas 6te choisi de caracteristiques provenant 
des arteres carotide's internes eloignees ou de l'artfere caro- 
tide primitive haute (qui serai t sujette a erreur de classi- 

10 fication, si elle etait s6lectionne<|, les seules donnees 

que l'on a reellement besoin de recueillir de chaque patient 
peuvent etre reduites aux donnees de variation de frequence 
Doppler de l'artfere carotide primitive basse (2 cm environ 
au-dessus de la ligne de partage) et de l'artere carotide 

15 interne proche f au lieu des donnees provenant des quatre 
emplacements qui ont 6t6 utilises dans les experiences de 
la presente invention. Ainsi, les donnees peuvent etre com- 
primees encore davantage. 

S'il presente plusieurs avantages sur les autres sys- 

20 t6mes f le format de donnees comprimees utilise dans la pre- 
sente invention peut ne pas etre le meilleur. Ce syst&me 
permet : (1) une evaluation sensible des variations rapides 
qui se produisent dans le signal Doppler avec le temps ; 
(2) une interpretation aisee du spectre ; (3) une represent 

25 tation adequate de l'asymetrie habituelle de la density spec- 
trale du signal ; (4) une fa9on sure de suivre le signal en 
presence de parasites en bande large et en bande etroite, 
d' intensity variable ; (5) une reduction de donnees suffi- 
sante pour faciliter la mise en memoire, le traitement et 

30 1* analyse ; et (6) des -donnees adequates pour 1* extraction 
et 1 • identification de caracteristiques instructives neces- 
saires a une bonne dif f Srenciation de l'etat d'avancement 
de la maladie. 

Une forme d'onde type pour l'ecoulement sanguin dans 

35 l'art&re carotide primitive proche d'une artere carotide 
interne normale se caracterise par un pic systolique pro- 
nonce , pointu et etroit, avec peu ou pas d'encoches dans 
la phase de deceleration , et une circulation diastolique 
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4 pic prononce. Une forme d'onde type pour l'ecoulement san- 
guin dans l'artere carotide primitive proche d'une artere 
carotide interne tres peu atteinte revele un pic systolique 
beaucoup plus bas, plus large et "arrondi", une encoche mar- 
quee au debut de la phase de deceleration, et une circulation 
diastolique plus amortie. Ces differences generales peuvent 
etre analysees categoriellement, evaluees et classees avec 
le format de donnees coraprimees et de reconnaissance de formes 
selon la presente invention. La premiere decision (normal 
ou malade) et la deuxieme (stenose superieure ou inferieure 
a 50 %) semblent importantes pour prendre une decision sur 
le traitement du patient. Mais les raisons qu'il y a a definir 
d'autres categories de stenose sont moins evidentes. On point 
de decision au premier niveau (15 ou 20 % d* occlusion) etait 
initialement base sur : (1) la necessite de graduer la maladie 
mieux que par une stenose simplement superieure ou inferieure 
a 50 %, pour permettre une evaluation de la progression du 
mal ; (2) le fait evident que 1 'accord entre radiologues 
lisant les mimes arteriographies est meilleur au niveau de 
20 % qu'au niveau de 10 % ; et (3) le fait qu'un minimum 
naturel se presentait dans la repartition de la base de don- 
nees de formation au niveau de 20 %. Accroitre la taille 
de la base de donnees peut permettre une meilleure gradation 
dans la categorie de 51 a 99 %. Bien qu'une meilleure discri- 
mination soit possible dans la categorie de 51 a 99 %, il 
ne semble pas realiste d'af finer encore le programme de faffon 
a distinguer des variations de 15 % environ en plus ou en 
moins, compte tenu du manque de precision de 1 • interpretation 
des arteriogrammes reels. 

Dans les experiences reelles, on s'est efforcd d'^valuer 
1' importance du technicien pour .1 'obtention d'un diagnostic 
correct quant au degre de stenose. Les essais de la presente 
invention consistent a prelever des donnees Doppler a ultrasons 
sur le meme patient dans trois h6pitaux diff brents, et k 
analyser ces donnees [obtenues par des techniciens differentsj 
avec l'ordinateur numerique dans le systeme de reconnaissance 
de formes qui est decrit dans cette invention. On a constate 
que dans environ 90 % descas, les diagnostics par ordinateur 
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sont les memes ou se retrouvent dans la meme categorie. Autre- 
ment dit, si nous classons le degre de stenose dans la cate- 
gorie de 36 a 50 % d 1 occlusion, dans 90 % des cas un techni- 
cien d'un autre hopital participant a l'essai trouvera que 

5 le degre de stenose est de 15 a 35 %, de 36 h 50 % ou de 
plus de 50 % (a l'interieur d'une meme categorie). Ainsi, 
le systeme s'est avere concluant en tant qu'outil non enva- 
hissant pour identifier la maladie et pour, tout au moins 
grossierement, specifier le degr£ de stenose. 

10 Le systeme de classification seloh la pr£sente invention 

est discret et ne tient pas compte des disaccords eventuels 
rencontres lors de chaque decision. II impose une decision 
pour une categorie. Par exemple, une artire class^e comme 
normale dans la premiere decision pourrait trfes bien, dans 

15 la seconde decision, etre classee comme prSsentant une ste- 
nose de plus de 50 %. Selon 1' ensemble de formation rSel 
et 1* ensemble d'essai envisage, il n'y a pas plus de 2,5 % 
de chances pour qu'une artfere presente plus de 20 % d 1 occlu- 
sion ou de stenose apres avoir et6 qualifi£e de normale dans 

20 la premifere decision. De meme, aucun vaisseau pr^sentant 
moins de 20 % de stenose n'a et6 classS comme ayant plus 
de 50 % de stenose avec le classif icateur r£el. Le classifi- 
cateur ne se trompe g6n£ralement que d'une categorie. 

Une erreur de classification r6sulte le plus souvent 

25 de stenoses tres prononcees chez des patients dont les donnees 
etaient caract^risees par un type de diagramme d'ecoulement 
quasi -permanent et non pulsatile. Ce type de diagramme d'^cou- 
leroent est mal reprSsente dans la base de donn£es de formation. 
A mesure qu'un plus grand nombre de donnees de ce type devien- 

30 dront disponibles, il ne fait pas de doute que l'algorithme 
de classification sera modifie, etant donn€ que la base de 
donnees de formation aura une meilleure base sur laquelle 
prendre les decisions de classification. Des diagrammes d'6cou- 
lement de ce type seront ajoutes h la base de donnees pour 

35 ameliorer celle-ci ainsi que la prise de decisions. 

Le niveau d' accord entre le diagnostic par ordinateur 
et 1» angiographic, tel qu'il est mesure par la statistique 
. Kappa, est egal a 0,77 ± 0,04 (erreur type). Ce niveau d'accard 
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se compare f avorablement avec le niveau d 1 accord qui peut 
etre atteint par des radiologues diff^rents lisant les memes 
art6riogrammes, aussi le syst&me apparait-il comme tres. pro- 
inetteur. 

Pour encore accentuer les distinctions et af finer le 
systeme a l'interieur de la categorie 51-99 %, une surveil- 
lance plus adequate de la progression de 1 'arteriopathie 
sera necessaire au cours d' Etudes de suivi & long terme. 
Pour determiner avec precision 1' emplacement d'^chantillon- 
nage £ l*interieur du vaisseau et pour r^duire la variability 
eventuelle impliguee par l'endroit o& se trouve le volume 
d'^chantillon Doppler a impulsions, on explorera £galement 
la possibility d'utiliser un dispositif de rep^rage des posi- 
tions tridimensionnelles. On s' attend & ce que cette m^thode, 
en permettant l'£chantillonnage de l'ecoulement a plus ou ■ 
moins 1 mm d'un emplacement deja echantillonne, trouve une 
application particuliere quand des etudes Duplex en s£rie 
sont effectuees. 

Le logiciel suivant represente le codage pr£f ere actuel- 
lement utilise pour effectuer les calculs necessaires pour 
analyser les caracteristiques au niveau des emplacements 
des art&res carotides . primitive et interne, et pour identifier 
des reperes important s dans les donnees, tels que systole, 
minima relatifs, ou maxima relatifs. Les observations rela- 
tives au logiciel et au codage FORTRAN classiques illustreront, 
pour l'homme du metier, le format de calcul utilise pour 
situer les caracteristiques. 
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35 



. TITRE MAXAMR 
. ENABL LC 
. GLOBL MAXAMR 

. GLOBL SRKMBA, PRCENT, MEDIAN, SWABTS, GTBYT 
. GLOBL SRTMIF 
MAXAMR — 

Sous-programme FORTRAN pouvant etre appele pour fournir 
une estimation resistante du mode (tend a ignorer les 
valeiirs ecarte*es, d'une fagon non lineaire). II fonc- 
tionne comiae suit : 

1) Trier le spectre sur les donnees d f entree, apres 
avoir mis l'ordre original dans 1' octet superieur de 
chaque mot. 

2) Prendre les elements "NROBUS" superieurs ou bien 
jusqu'a ce que le niveau "BRUIT" soit atteint, ou bien 
jusqu'a ce que le pourcentage "LIMITE" du maximum global 
soit atteint r et 

3 ) . . • Fournir la m^diane de ce groupe en tant que mode . 
A 1" entree — 

Le fichier d" entree est suppose Sfcre positionne au debut 
du spectre (ou autre) a examiner, 

2(R5) est I'adresse de la rangee comportant le spectre. 
64 (R5) est le nombre d 1 elements (mots) dans la rangee. 
66 (R() est le H degr3 de robustesse" voulu. II d^terminera 
le nombre maximal de points de donnees a considerer pour 
effectuer 1 'estimation. 
A la sortie — 

L'estimation en nombre entier est retourn£e (58. (R5) . 
Le fichier d" entree est avance a l 1 octet suivant la 
rang£e qui vient d'Stre analysee. 
Constantes d'assembleur. 
LIMITE =50. 
BRUIT =5. 
MAXAMR : 

; Forcer a des valeurs raisonnables pour les parametres 
; d' entree. 

TST 64 (R5) 

BGE l<j 
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1«: 



CLR 

TST 
BGE 
CLR 



fc4(R5) 
06(R5) 
C6(R5) 



2*: 

; Forcer la valeur "R M i R4, o6 elle restera. 

MOV ©6(R5), R4 

; Mettre des index de frequence dans les octets superieurs 
de la rangee d f entree. 

MOV 2<R5), Rl ; Rlj = adresse de la rangee 

MOV ©4CR5), R2 ; R2: . nombre d'Sl&nents de la 

rangee. 

CLR R3 ; R3: « index de frequence. 



SETKY 



MOVB 

SWAB 

BIS 

SWAB 

INC 

MOV 

SOB 

MOV 



(Rl),RO 
RO 

R3,RO 

RO 

R3 

RO, (Rl)+ 
R2 f SETKY 
ML, "(SP) 



; Trier la rangSe du plus petit au plus grand. 

MOV 2(R5),R1 

MOV <S 4 (R5),R2 

APPELER SRKMBA 
; Calculer 1' amplitude de troncature 



donn£e suivante & RO 

mettre octet sup6rieur en position 

masque dans 1' index de frequence 

retablir ordre 

frequence suivante 

mettre en mdmoire rdsultat dans 

rangee 

tourner sur boucle jusqu'S ce 
que ce soit fait 
sauvegarder le point final. 



MOV 
DEC 
MOVB 
MOV 

MOV 

APPELER 



(SP),R2 
R2 

-ll.(R2),Rl 

# BRUIT, R3 

# LIMITE , RO 
PRCENT 



Retablir fin de rangee 



Rl: 
R3: 



amplitude maximale 
niveau au-dessous duquel 
les donnees sont du bruit 



Descendre en suivant la liste. 



RO: = amplitude de seuil. 
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; R2: amp, max./8 
; augmenter M R" de (amp. max.)/8 
; R2: = nombre de points a considSrer 
? depart au maximum 



; au-dessous du bruit ? 
; oui 

; au-dessous du seuil ? 
; oui 



MOV -1(R2),R1 

MOVB Rl,R2 

ASH #-2,R2 

ADD R2 ,R4 

MOV R4,R2 

MOV (SP)+,R1 
DESCND : 

TST -(Rl) 

CMPB (R1),R3 

BLE ESTIMA 

CMPB (Rl), ro 

BLE ESTIMA 

SOB R2, DESCND 

; Calculer 1 'estimation. Rl adresse le point oil la descente 
est arretee. R4 +1 -R2 donne le nombre de points dans 1» esti- 
mation. 

SUB R4,R2 

NEG R2 ; R2: = nombre de points dans 

1 1 estimation 

INC R2 ' s f assurer qu'il y en a au moins un 

APPELER SWABTS ; permuter toutes les paires d' octets 

Trier en fonction de frequence 

APPELER SRKMBA 
Fournir la mSdiane 

APPELER MEDIAN 

BIC ^177400^0 

MOVB RO,RO 

MOV RO,68. (R5) 

RETOUR 

. PIN 



35 



2591884 



. TITRE LIMITES 

. ENABL LC 

• GLOBL LIMITES 

. GLOBL SRKMBA, MEDIAN, SVABTS, WINEST 
; LIMITES ~ 

; Calcule une estimation "robuste" des contours haut et 
; bas d'un spectre autour du point d f implantation. 
; A l f entree — 

; G2(R5) = Estimateur d 9 emplacement 

(R5) = Niveau de seuil ou de contour 
; 6(R5) = Adresse la" rang£e nombre entier avec le spectre 
; ayant des index de frequence dans les octets 

superieurs de chaque Element. 
; 68 (R5) = Nombre d 1 el&nents dans la rangee du spectre. 
; @io(R5) = Robustesse. 
; A la sortie — 

Cl2(R5) = limite inferieure 
? Gl4(R() = limite sup^rieure 
; Constantes d'asserableur . 
BRUIT = 1 
LIMITES : 

; Mettre le seuil dans R3 
MOV G4(R5),R3 

CMP R3,$ BRUIT ; seuil dans marge de bruit ? 

BGE 1$ ; non 

MOV # BRUIT, R3 ; oui, hausser le seuil. 

; Calculer adresse du point d' implantation. 
MOVB 62(R5),R1 

ASL Rl ; transformer en nombre d 1 octets. 

ADD 6(R5),R1 

MOV R1,R2 ; Rl et R2 adressent tous deux 

le point d ' implantation 
; Rechercher & partir du point d ' implantation en descendant 
; pour trouver la limite inferieure 

MOVB G2(R5) / R4 

BEG LOWDNE 

TST (Rl)+ 
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ADD 
LOWSEA: 



42, R4 



DEC 


R4 


BLT 


LOWDNE 


DEC 


Rl 


CMPB 


-(Rl),R3 


BGE 


LOWSEA 



R4 = Distance a moins I'infini. 



rechercher la fin du tampon 



LOWFND: 

; Controler la coherence dans le point suivant. 



10 



MOV . 
DEC 
CMPB 
BGE 
LOWDBE: 



R1,R0 
R0 

-(R0) r R3 
LOWSEA 



point suivant < seuil ? 
non, continuer a chercher. 



15 



? Rechercher k partir du point d 1 implantation en montant 
; pour trouver la limite sup^rieure 





MOV 


&8. (R5),R4 








MOVB 


G2(R5),R0 








SUB 


RO,R4 






20 


SUB 


#2,R4 




R4: = Distance a plus I'infini. 




BEG 


HIDNE 








TST 


-(R2) 








HISEA: 










DEC 


R4 






25 


BLT 


HIDNE 


• 


rechercher la fin du tampon. 




TST 


(R2) + 


; 


avancer 1' index au mot suivant. 




CMPB 


(R2),R3 


; 


au-dessous du seuil ? 




BGE 


HISEA 


« 


non, contr61er le mot suivant. 




HIFND: 








30 


Controler la coherence dans le point suivant. 




DEC 


R4 








BLT 


HIDNE 


r 


rechercher la fin du tampon. 




TST 


(R2) + 


• 
r 


avancer 1" index au mot suivant 




CMPB 


Gr2,R3 


t 


point suivant < seuil ? 


35 


BGE 


HISEA 


7 


non, continuer a chercher. 




TST 


* -(R2) 


/ 


oui, corriger pour extra-increment. 



HIDNE: 

; Transformer Rl et R2 en valeurs d' index. 
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SUB 6(R5),R2 

SUB 6{R5),R1 
Fournir limites inf erieure et superieure. 

ASR R1 ; transformer d' octet en. 

ASR R2 .... nombre d' Elements. 

MOV Rl,6l2.(R5) 

MOV R2,ei4. (R5) 

RETOUR 

.PIN 
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SOUS -PROGRAMME "REPERES" (BUFF, NPTS, NDELX, NX, NY, 
XSCALE, YSCALE , et NlMAXl, NlMINl, N1ZR01, RlARBl, NSYSTO, 
R1ARE2 ) 

OCTET BUFF (1) 

5 c Calcule une certaine information de repere int£ressante 
C h partir de la sequence "BUFF" qui est suppos£e contenir 
C une premi&re fonction dSrivee. 

C Trouver le point d 'acceleration maximale dans la region 
systolique. 

10 APPELER ' EXTB YT ( BUFF (50), 120, NlMINl, NlMAXl) 

NlMINl - N1MIN+50 
N1MAS1 = NlMAXl+50 
C Trouver le pic systolique suppose. Celui-ci est defini comrae 
C le premier point de pente proche de zero apr&s NlMAXl. 
15 FAIRE 400, NSYST0=N1MAX1 , NlMINl 

SI ( (BUFF(NSYSTO). LE. BUFF ( NlMINl )+l) 
ft .OU. ( BUFF ( NS YSTO ) .LT. 1) ) ALLER A 500 
400 CONTINUER 
500 CONTINUER 

20 C Trouver d'autres extremes maintenant que NSYSTO est 
dSfini 

NSEARC - NSYSTO+10 

APPELER EXTBYT (BUFF (NSEARC) , 6 5, NlMINl, NlMAXl) 
N1MIN1= N1MIN1+NSEARC 
25 N1MAX1= N1MAX1+NSEARC 

APPELER EXTB YT( BUFF ( (NPTS/2 )+l ) , NPTS/2, N1MAX2) 
NlMIN2= N1MIN2+NPTS/2 
N1MAX2= N1MAX2+NPTS/2 
C Trouver le premier point de pente zero apr&s le point 
30 C de deceleration maximale (NlMINl), que l'on espfere etre 
C le premier point de pente zero aprfes la systole. 
FAIRE 600, N1ZR01=N1MIN1, 200 

SI (BUFF(NIZROI) .GT. -1) ALLER A 700 
600 CONTINUER 
35 700 CONTINUER 

C Calculer l'aire sous la courbe de deceleration post- 
C systolique au point de deceleration maximale. 
R1ARE1=0 



39 



2591884 



R1ARE2=0 

FAIRE 300, I=N1MIN1-10, NlMINl+10 

II=BUFF (I) 

SI (II .LT. 0) RlARE2= R1ARE2 + II 
5 300 R1ARE1=R1ARE1+BUFF(I) 

R1ARE1= Rl ARE1/FLOAT ( NDELX ) 

R1ARE2= Rl ARE2/PLOAT ( NDELX ) 
C Bcrire barres sur la courbe pour que les gens puissent 

voir ce qui s'est pass6. 
10 APPELER BRAKIT ( NX , NY , NlMINl , NlMAXl ,90, SXCALE , YSCALE ) 

APPELER BRAKIT (NX, NY, NS YSTO , Nl ZROl ,100, XSCALZ , YSCALE ) 

FIN 
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SOUS- PROGRAMME BFEATU( LOCBUF, SCLBUF, 

& NPTS, YSCALE , NXO , NYO , CBUFF 1 , CBUFF2 , NROBUS , NCYCMX,PA3NAM) 

COMMUN/REPERES/NOMINl , NOMAXl , NOMIN2 , NOMAX2 , NlMINl , NlMAXl , 

& NlMIN2,NlMAX2,NSYSTO,NlZR01 

COMMU N / P ARAMS /A IRE , COSANG ( 3 ) , SINANG { 3 ) , 

& R0VAL1 , R0VAL2 , R0ARE1 ,R0ARE2 , 

& Rl VALl , Rl VAL2 , Rl ARE1 , R1ARE2 , 

& RMVAL1 , RMVAL2 , RMAREl , RMARE2', 

& RPVALl , RPVAL2 , RPARE1 , RPARE2 , 

& PENTE 

COMMUN/RANGEES/SCR1 , SCR2 ,SCR3, RTIME , RMAG , RPHASE 
CBUFF1(1),CBUFF2(1) VIRTUELS 
CBUFF1, CBUFF2 , CTEMP COMPLEXES 

RTIME ( 512 ) , RMAG ( 512 ) , RPHASE ( 512 ) , FEATS ( 1 ) REELS 
OCTET PATNAM ( 1 ) 

OCTET LOCBUF ( 1 ) , SCLBUF ( NPTS , 1 ) 

OCTET SCRl ( 512 ) , SCR2 ( 512 ) r SCR 3 ( 512 ) , NOTCHS ( 2 ) 

PI2= 3,14159265*2 

NPTP2=512 

XECHELLE= YECHELLE/ 3 . 

Copier les donnees d 1 emplacement pour ef facer les tam- 
pons 1, reels, et complexes. 
APPELER ST COPIE ( LOCBUF , SCRl , NPTS ) 

Forcer la fin reelle des donnees a se produire a t=0,5 
sec. 

APPELER INIBYT(LOCBUF(240),41,LOCBUF(240) 
APPELER INIBYT( SCLBUF (240,1) , 41, SCLBUF (241,1) ) 
APPELER INIBYT ( SCLBUF { 24 0 , 2 ) , 41 , SCLBUF (240,2)) 
APPELER INIBYT ( SCLBUF ( 2 4 0 , 3 ) , 4 1 , SCLBUF (240,3)) 
APPELER INIBYT ( SCLBUT ( 2 4 0 , 4 ) , 4 1 , SCLBUF (240,4)) 
Soustraire le decalage DC dans les donnees d' entree. 
APPELER ADDBYT( SCRl, NPTS, -30) 
APPELER ADDBYT( LOCBUF, NPTS, -30) 
FA I RE 100,1=1,4 

APPELER ADDBYT ( SCLBUF (1,1), NPTS , -30 ) 
Enlever la pointe de declenchement de l'ECG 
SCRl (41)= SCRl (40) 
LOCBUF (41)= LOCBUF (40) 



41 



2591884 



FAIRE 110, 1=1,4 

110 SCLBUF (41,1)= SCLBUF(40,I) 

APPELER DOTIME ( LOCBUF , SCLBUF , NPTS , RTIME , NPTP2 , AIRE; 
PATNAM) 

SI (AIRE .GT. O.O,) ALLER A 120 
IMPRIMER 111, (PATNAM(I),I=1,8) 

111 FORMAT ( 1 BFEATU-W-',8A1, "a <*nergie nulle") 
ALLER A 5000 

120 CONTINUER 

C Tracer la courbe des donn^es d 1 emplacement . 

APPELER PLTRE1 ( NXO r NYO , O , 1 FIX (30 * YECHELLE ) , 

& XECHELLE , YECHELLE , RTEMPS ,280 / NROBUS , NC YCMX , PATNAM ) 
C Initialiser les coordonn^es de tracage. 

NXLO=NXO+NPTS+4 0 

NYLO=NYO- 30*YECHELLE 

APPELER SCAIRE ( RTEMPS ,512,1. /AIRE ) 

APPELER REACOM ( RTEMPS , NPTP2 , CBDFF1 ) 
C Trouver des repfcres dans la forme d'onde de Vitesse. 
C (Otiliser les donnees "SCR1" pour les extr&nes Stant 

donn<5 qu'elles ont deS donnees de + 6 sec. intactes). 

APPELER EXTB YT ( SCR1 , (NPTS/2 )+30,NOMINl,NOMAXl ) 

APPELER EXTBYT(SCR1( <NPTS/2)+l ) ,NPTS/2,NOMIN2,NOMAX2) 

NOMIN2= NOMIN2+ NPTS/2 

NOMAX2= NOMAX2+ NPTS/2 

APPELER BRAKIT(NXO,NYO,(NOMAX1-1), (NOMIN2-1 ) ,35, XECHELLE, 
YECHELLE) 

: Calculer la premiere ddriv^e et laisser dans le brouillbn 3. 
NDELX-9 I Instaurer delta X. 
APPELER DERI V ( LOCBUF , NPTS / NDELX , SCR2 ) 
APPELER LISSE(SCR2,NPTS,4,SCR3) 

APPELER REPERES ( SCR3 , NPTS , NDELX, NXLO , NYLO, 1 . , YECHELLE, 

& N1MAX1 , N1MIN1 , Nl ZROl , R1ARE1 , NS YSTO , R1ARE2 ) 

APPELER BRAKIT ( NXO , NYO , NS YSTO , Nl ZROl ,40, XECHELLE , YECHELLE ) 

APPELER SCABYT(SCR3,NPTS,2, 5) 

APPELER ADDBYT(SCR3,NPTS,50) . 

APPELER PLTBYT(SCR3,NPTS,NXLO,NYLO,l. , YECHELLE) 
NXFIL— NXO 

NYFIL= NYO- 1 0 Q* YECHELLE- 2 0 
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C Calculer la DFT et afficher la grandeur. 
APPELER HAMEND { CBUFF1 , NPTP2 ) 
APPELER FFTFOR ( NPTS , CBUFF1 ) 
C Mettre en phase sur "RPHASE" 
5 FA IRE 200, I=1,NPTP2 

200 RPHASE(I)= AIMAG( CLOG( 1.0E-6+CBUFFl(I) ) ) 

C Annoncer la phase. 
FAIRE 400, 1=1,511 

NDIFF= (ABS(RPHASE{I+1)-RPHASE(I))/PI2 ) +0,5 
10 IF ( RPHASE (1+1) .LT. RPHASE(I) ) NDIFF= -NDIFF 

400 RPHASE ( I +1)=RPHASE(I+1) - PI2*NDIFF 

PENTE= -2*RPHASE( (NPTP2/2)+l ) 
FAIRE 450, 1=1, NPTP2 
450 RPHASE (I )= RPHASE(I) + PENTE* ( 1-1 ) /NPTP2 

15 C APPELER REAREA ( RPHASE , NPTP2 , RMA6 ) 

C APPELER EXTREA ( RPHASE , NPTP2 , RMIN , RMAX ) 

C APPELER SCAREA(RPHASE,NPTP2,50./(AMAX1(ABS(RMAX),ABS 

C (RMIN))) ) 

APPELER ADDREA ( RPHASE ,NPTP2 , 50 . ) 
20 C NXLO=NXLO+NPTS+40 

C APPELER PLTREA ( RPHASE ,NPTP2 , NXLO , NYLO 0,55, YECHELLE ) 
C APPELER REAREA ( RMAG , NPTP2 , RPHASE ) 
C Et maintenant pour la grandeur. 

APPELER COMMAG ( CBUFF1 , NPTS/4 , RMAG ) 
25 APPELER AREREA ( RMAG ( 6+1 ) ,5,RMARE1) 

APPELER EXTREA ( RMAG , NPTS/4 , RMIN , RMAX ) 
APPELER SCAREA( RMAG, NPTS/4, 100. /RMAX) 
NXLO= NXLO+NPTS +40 
NYLO= NYLO - YE CHELLE * 1 0 0 - 30 
30 APPELER BRAKIT{NXLO,NYLO,5,9,30,8., YECHELLE 

APPELER PLTREA ( RMAG , NPTS/8 , NXLO , NYLO , 8 . , YECHELLE 
C Laisser l^rigine de la courbe en has a gauche de la 

rangee suivante. 
5000 NXO= 30 
35 NYO= NYO- YECHELLE * 1 0 0 - 60 

D TYPE *, l COSANG[l,2,3J = » , (COSANG( I) , 1=1,3) 
D TYPE *, 'NSYSTO= f ,NSYSTO 
FIN 
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SOUS -PROGRAMME BMAKC1 ( FEATUR , LOCBUF , SCLBUP , NPTS , I FEAT/ 
NMAX , ITER AT , DIALOG ) 

c edcece^eeeeecec&tccecc^c teeter t&ct 

C D^finit et calcule toutes les caract^ristiques pour 

5 C 1 * emplacement n°l de carotide primitive. 

C A 1' entree — 

C COMMUN hloque DATAP,REPERE§PARAMS, et RANGEES contient 

C des informations valables pour cet emplacement. On peut 

C les ddfinir par un appel au sous-programme BFEATU . 

10 C LOCBUF, SCLBUF ont 1 ' information d 1 emplacement et d'^chelle 
pour cet emplacement. 

C NPTS= nombre de points valides dans LOCBUF/SCLBUF 

C FEATURE rangee reelle dans laquelle les caract£ristiques 

C ("features") seront ecrites. 

15 c IFEAT= d£calage dans la rangee FEATUR pour la premiere 
caract^ristique de ce module. 

C NMAX= dimension de la rangee FEATUR pour controler les 

C bornes de la rangee. 

C A la sortie — 

20 C SI il reste assez de place dans FEATUR 

C les nouvelles caract£ristiques de ce module y auront 

C ete laiss^es. 

C IFEAT= IFEATl & 1' entree ]+ (nombre de nouvelles 

caractdristiques introduites ) . 

25 c SINON, sortir. ^ ^, * „ & * a & 

RTIME ( 512 ) , RMAG ( 512 ) , RPHASE ( 512 ) , FEATUR ( 1 ) REELS 
OCTET LOCBUT ( NPTS ) , SCLBUF ( NPTS , 1 ) , SCR1 ( 512 ) , SCR2 ( 512 ) , 
SCR3(512) 

30 DIALOGUE LOGIQUE, PUPITRE . 

COMMUN/DATAP/NVERSI , NVRMAX , NVRSCL , NCYCMX / NCYCSC 
COMMUN/REPERES/NQMIN1 , NOMAXl , NOMIN2 , NOMAX2 r N1MIN1 , N1MAX1 , 
& N1MIN2,N1MAX2,NSYST0,N1ZR01 

COMMUN / INDMRK/NAS YS T , NAlZRl , NBSYST , NB1ZR1 , NCSYST , NCI ZR1 
3 5 COMMUN/PARAMS/AIRE , COSANG ( 3 ) , SINANG ( 3 ) , 

& ROVALl r R0VAL2 , ROAREl , ROARE2 , 
& Rl VAL1 , Rl VAL2 , Rl AREl , R1ARE2 , 
& RMVAL1 , RMVAL2 e RMAREl , RMARE2 , 
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& RPVALl / RPVAL2 , RPARE1 , RPARE2 / 
& PENTE 

COMMUN/RANGEES/SCRl , SCR2 , SCR3 , RTIME, RMAG , RPHASE 
COMMUN/BMACK1/AMIMOD , AMXMOD , AMI UP 8 , AMXYP8 , POSTDC 

C Constantes de module 

ELON= l.OE-10 ! Zero de machine 

NUMFEA= 50 i Nombre de caracteristiques dans 

BMAKC1 

10 IBR01=22 J Index de la premi&re caract£ris- 

tique d' elargis semen t spectral. 

SI ( ( ITERAT .EQ. O) .ET. DIALOG) . 

& IMPRIMER *, 'BMACKl-I- Decalage d' index de caracteris- 
15 tiques = * , IFEAT 

C S 1 assurer qu'il y a asse2 de place dans le tampon. 
SI { ( IFEAT+NUMFEA) .LE. NMAX) ALLER A 10 
SI (ITERAT .EQ. O) 

& IMPRIMER *, 'BMAKCl-F-Tampon trop petit. 1 
20 RETOUR 

10 CONTINUER 

APPELER INIREA ( FEATUR , NUMFEA# O.O) 
SI (AIRE .EG. O.O) ALLER A 5000 
C Calculer les caracteristiques de la forme d'onde de 
25 vitesse. 

C Etablir des rep&res pour cette art&re. 
NASYST= NSYSTO 
NA1ZR1= NIZROl 

C DDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDD 
30 C Ddfinir toutes les caracteristiques de cette art&re 

le mieux possible. 

FEATUR { 1 +1 FEAT ) » AIRE 

FEATUR ( 2 +1 FEAT ) = COSANG(l) 

FEATUR ( 3+IFEAT ) - SINANG(l) 
35 C Donnees de contour a la systole. 

C Aire sous le pic systolique (+/- 15 intervalles) • 

SYSARE= FFARER ( RTIME , NPTS , NSYSTO- 15 / NSYSTO+15 ) 

FEATUR (4+IFEAT)= SYSARE*AREA/COSANG ( 1 ) 
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Valeurs de mode. 

FEATUR ( 5+IFEAT)= FFAVGB ( LOCBUF , NPTS ,20,5) i Fin de diastole 
FEATUR( 6+1 FEAT) = LOCBUF { NSYSTO ) ! Pic systolique. 
FEATUR ( 7+IFEAT ) = LOCBUF ( Nl ZROl ) i Premier point de pente 

z6ro 

FEATUR ( 8 +IFEAT ) = LOCBUF ( NOMIN2 ) i Encoche dicrotique. 
Degr6 de deceleration suivant la systole. 
FEATUR ( 9 +IFEAT ) =R1ARE1/C0SANG ( 1 ) 
Ddbit diastolique 

FEATUR (10+1 FEAT )= RTIME IMIKO { NOMIN 2 +50,280)) *AIRE COSANG(l) 
FEATUR ( 11+IFEAT ) = RTIME (20) * AIRE/COSANG ( 1 ) 
Frequence d' amplitude maximale dans la WFT entre les 
harmoniques 4 et 9. 

APPELER EXTRE2(RMAG,NPTS,4,9, IMIN, IMAX> 
FEATUR (12+IFEAT)= IMAX 

Frequence d' amplitude maximale dans la WFT entre les 
harmoniques 5 et 9. 

APPELER EXTRE2(RMAG,NPTS,5,9,IMIN,IMAX1) 
FEATUR ( 1 3 +IFEAT ) = IMAX1 ! Caracteristique originale 13. 
Frequence d f amplitude maximale dans la WFT entre les 
harmoniques 6 et 9. 

APPELER EXTRE2(RMAG,NPTS,6,9,IMIN,IMAX) 
SI (IMAX .BQ. (IMAXl+1)) IMAX=IMAX1 

FEATUR (14+1 FEAT ) = IMAX 1 Version affile de l'ancienne 

caracteristique 13. 
Frequence d' amplitude minimale dans la WFT entre 2 et 
le maximum precedent. 
I=IMAX 

APPELER EXTRE2 ( RMAG , NPTS ,2,1, IMIN , IMAX ) 
FEATUR ( 1 5 +IFEAT ) = IMIN 

Index de pulsatilite, rapport magique, etc. 

POSTDI= FFAVGB ( SCLBUF (1,4), NPTS , 20 , 10 ) 

SI (POSTDI .LE. 0) P0STDI=» 1 

FMEAN=* FFAVGR(RTIME, NPTS, 160,120) 

SI (FMEAN. LE. O) FMEAN» 1 

FSYSTO= RTIME (NSYSTO) 

FEATUR (15+1 FEAT )= (FSYSTO-RTIME(NlZROl)/ 

& FMEAN i Rapport de Pourcelot 
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C Rapport d'Yve 

FEATUR ( 16+1 FEAT ) = ( FLOAT ( SCLBUF ( NSYSTO , 4 ) -SCLBUF 
& ( Nl ZROl , 4 ) ) / { SCLBUF ( NSYSTO , 4 ) - 2 9 , 5 ) ) + 

0,5 

FEATUR ( 1 7 +IFEAT ) = FSYSTO/POSTDI 
C Ce qui suit reflete la presence ou !• absence d'une 
"fenetre systolique" . 

APPELER EXTBYT (SCLBUF (1,1), NPTS , NS YSTO-1 2 ,20, IMIN , IMAX > 

FEATUR ( 1 8+IFEAT ) = SCLBUF ( IMAX, 1 ) -30 ! 3 dB inffirieur. 

APPELER EXTBYT ( SCLBUF (1,3), NPTS , NS YSTO-1 2 ,20, IMIN , IMAX ) 

FEATUR ( 1 9+IFEAT )= SCLBUF (IMAX, 3) -30 ! 9 dB infSrieur. 
C Information de phase. 

FEATUR ( 2 0+IFEAT ) = RPHASE ( 2 ) 

FEATUR ( 21+IFEAT ) = RPHASE ( 3 ) 

FEATUR (22+IFEAT)= RPHASE (4) 
C CaractSristiques bashes sur l'elargissement spectral 

(il y en a actuellement 28). 
C Fin de diastole. 

FAIRE 800, 1=1,4 

800 FEATUR(IBROl+0+I+IFEAT)= 

& FADIFB(LOCBUF,SCLBUF(1,I),NPTS,20,5) 
C Systole. 

900 PA FEAT^tl5Ril+4+I+IFEAT) = 

& FADIFB ( LOCBUF , SCLBUF (1,1), NPTS , NSYSTO , 5 ) 

C Premier point de pente zero. 

FAIRE 1000, 1=1,4 
1000 FEATUR ( IBR01+8+I+IFEAT) = 

& FADIFB ( LOCBUF , SCLBUF (1,1), NPTS , NlZROl , 5 ) 
C 100 msec, apr&s la systole. 

FAIRE 1600, 1=1,4 
1600 FEATUR ( IBR01+12+I+IFEAT ) = 

& FADIFB ( LOCBUF , SCLBUF (1,1), NPTS , NS YSTO+4 0,5) 

FAIRE 1650, 1=1,4 
1650 FEATUR ( IBR01+16+I+IFEAT ) = 

& FADIFB ( LOCBUF , SCLBUF (1,1), NPTS , NSYSTO , 5 ) /COSANG ( 1 ) 

FAIRE 1700, 1=1,4 
1700 FEATUR (IBROl+20+I+IFEAT)= 

& FADIFB ( LOCBUF, SCLBUF (1,1), NPTS , NS YSTO+30 ) /COSANG ( 1 ) 
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FA I RE 1900, 1=1,4 
1900 FEATUR(IBR01+24+I+IFEAT)= 

& FAD I FB ( LOCBUF , SCLBUF (1,1), NPTS , NS YSTO+4 0 ) /COSANG ( 1 ) 
C Enregistrer certaines valeurs irit£ressantes pour la 

posterity. 

P0STDC= POSTDI ! Fin d^coulement diastolique. 
C DDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDD 
C IFEAT= nombre de caract^ristiques actuelles, et est 

egalement laissd en Evidence a la fin de la liste des 

caracteristiques • 
5000 IFEAT= IFEAT+NUMFEA 

SI (NUMFEA .NE. (IBROl+28) 

& IMPRIMER *, "Les caracteristiques dSfinies et engendrees 

par BMACK1-W ne correspondent pas". 

FIN 
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SOUS -PROGRAMME BMAKI 1 ( FEATUR , LOCBUF , SCLBUF , NPTS , I FEAT , 
NMAX , ITER AT # DIALOG) 

C D^finit et calcule toutes les caracteristiques pour 
5 c !• emplacement n° 1 de la carotide interne. 
C A 1' entree — . 

C COMMUN bloque DATAP , REPERES , PARAMS , et RANGEES contient 
C des informations valides pour cet emplacement. On peut 
C les definir en appelant le sous-programme BFEATU. 
10 c LOCBUF, SCLBUF ont 1 1 information d 1 emplacement et d'Schelle 

pour cet emplacement. 
C NPTS= nombre de points valides dans LOCBUF/SCLBUF 
C FEATURE rangee r6elle dans laquelle les caracteristiques 

seront 6crites. 

15 c IFEAT= decalage dans la rangee FEATUR pour la premifere 
caracteristique de ce module. 
C NMAX= dimension de la rangee FEATUR pour controler les 

borne s de la rangee. 
C A la sortie — 
20 c SI il reste assez de place dans FEATUR 

C les nouvelles caracteristiques de ce module y auront 
C 6te laiss€es. 

C IFEAT= I FEAT 1 9 entree ] + (nombre de nouvelles caract£- 
C ristiques introduites) 
25 C SINON, sortir. 

RTIME ( 5 1 2 ) , RMAG (512), RPHASE ( 512 ) , FEATUR { 1 ) REELS 
OCTET LOCBUF ( NPTS ) r SCLBUF ( NPTS , 1 ) , SCR1 ( 512 ) , SCR 2 ( 512 ) , 
SCR3(512) 
30 DIALOGUE LOGIQUE, PUPITRE 

COMMUN/DATAP/NVERS I , NVRMAX , NVRSCL f NC YCMX , NC YCSC 
COMMUN/REPERES/NOMIN1 , NOMAX1 , NOMIN2 , NOMAX2 f NOMIN1 , N1MAX1 , 
& N1MIN2,N1MAX2,NSYSTO,N1ZR01 

COMMUN/INDMRK/NASYST , NA1 ZR1 , NBSYST , NBlZRl f NCSYST , NCI ZRl 
3 5 COMMUN/PARAMS/AIRE , COSANG ( 3 ) , SINANG ( 3 ) , 

& ROVAL1 , ROVAL2 , ROARE1 , ROARE2 , 
& R1VAL1,R1VAL2,R1ARE1,R1ARE2, 
& RMVAL1 , RMVAL2 , RMARE1 , RMARE2 , 
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& RPVALl , RPVAL2 , RPARE1 , RPARB2 , 
& PENTE 

C0MMUN/RANGEES/SCR1 , SCR2 , SCR3 , RTIME , RMAG , RPHASE 
COMMUN/ AMAKC1 /AMIMOD , AMXMOD , AMIUP8 , AMXOP8 , POSTDC 
5 COMMUN/BMAKIl/BMIMOD,BMXMOD,BMIOP8,BMXOP 8,MAX3DI,MAX9DI 

C Constantes de module. 
ELON=* 1. OE-10 * 
NUMFEA= 44 
10 IBR01= 16 

SI ( ( ITERAT .EQ. 0) .ET. DIALOG) 

IMPRIMER *, BMAKIl-I-Decalage d' index de caracteris- 
tiques = I PEAT 

15 c S 1 assurer qu'il y a assez de place pour le tampon. 
SI ( (IFEAT+NUMFEA) .LE. UMAX) ALLER A 10 
SI (ITERAT .EQ. O) 

& IMPRIMER *, BMAKIl-F-Tampon trop petit. 

RETOUR 
20 10 CONTINUER 

C FAIRE 50 , 1=1 , NPTS 
C50 Type 51,I,LOCBUF(I) , - 

C & SCLBUP(I f l) f SCLBOF(I,2),SCLBUF(I # 3) f SCLBUF(I r 4) 
C51 FORMAT ( 1 ' , 13 ,X, 5 (14 r X) ) 
25 C APPELER INIREA ( PEATUR { IFEAT+1 ) , NUMFEA, O.O) 
. SI (AIRE .EQ. O.O) ALLER A 5000 
C Calculer les caracteristiques de la forme d'onde de 
Vitesse. 

C Etablir des reperes pour cette artere. 
30 NBSYST= NSYSTO 

NB12R1= NlZROl 

C DDODODDDDDDDDDDDDDODDDDDDODDDDDDDDODDDDDODODDDODDDDOODDD 
C Definir toutes les caracteristiques* 
FEATUR(1+IPEAT)= AIRE 
35 FEATUR ( 2+IFEAT ) = C0SANG(2) 

ADP= FFAVGB(SCLBUF<1,4), NPTS, 10,10) ! Frequence de 

crete de fin de diastole. 
SI (ADP .LE. 25) ADP= 0 i Imposer 0 pour < = 2,5 kHz. 
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FEATUR ( 3+IFEAT ) = ADP 

SI (ADP .LE. 40) ADP= 0 ! Imposer 0 pour < = 4 kHz. 
C Aire sous le pic systoligue (+/- 15 intervalles ) . 
C SYSARE= FFARER(RTIME,NPTS,NSYSTO-15,NSYSTO+15) 
5 C FEATUR ( 4 +IFEAT)= SYS ARE* AI RE/COS ANG { 2 ) 
C Frequence de crete de fin de diastole. 

FEATUR { 5+IFEAT ) = FFAVGB ( SCLBUF ( 1 , 4 ) ,NPTS, 10,10) ! Pin de 

diastole. 

FEATUR (6+IFEAT)= R1ARE2/COSANG ( 2 ) 
10 c Frequences maxima les. 

APPELER EXTBYT ( SCLBUF ( 1 , 2 ) , 1 6 0 , IMIN , IMAX ) 
MAX3DI= SCLBUF ( IMAX, 2 ) 

FEATUR ( 7 +IFEAT ) =*MAX3DI ! Frequence 3 dB de crSte 

APPELER EXTBYT ( SCLBUF (1,4), 160 , IMIN, IMAX) 
!5 MAX9DI= SCLBUF (IMAX, 4) 

FEATUR (8+IFEAT)= MAX9DI ! Frequence 9 dB de crete 

C Degr6 de deceleration suivant la systole. 

FEATUR ( 9+IFEAT ) = R1ARB1/COSANG { 2 ) 
C Energie dans I 1 amplitude raaximale entre les harmoniques 
20 5 et 9. 

APPELER EXTRE2 ( RMAG , NPTS ,5,9, IMIN , IMAX ) 
FEATUR ( 1 0+IFEAT j = IMAX 
C Frequence d' amplitude minimale dans la WFT de 2 au maxi- 
mum precedent. 
25 I«IMAX 

APPELER EXTRE2(RMAG,NPTS,2,I,IMIN,IMAX) 
FEATUR (11+IFEAT)= IMIN 
C Frequence d 1 amplitude maximale entre les harmoniques 
4 et 9. 

30 APPELER EXTRE2(RMAG,NPTS,4,9,IMIN,IMAX) 

C Energie totale dans les harmoniques 8 k 12. 

ENERGIE=FFARER ( RMAG , NPTS ,8,12) 
C Index de pulsatilite, rapport magique, etc. 
POSTDI=* FFAVGB ( SCLBUF (1,4), NPTS ,20,10) 
35 SI (POSTDI .LB. 0) POSTDI=l 

FMEAN=* FFAVGR(RTIME, NPTS, 160, 120) 
SI (FMEAN .LE. 0) FMEAN= 1 
FSYSTO=* RTIME ( NS YSTO ) 
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FEATUR(12+IFEAT)= (FSYST0-RTIME(NlZR01 ) )/ 

& FMEAN ! Rapport de Pourcelot. 

C Rapport d'Yve. 

FEATUR (13+1 FEAT ) = ( ( SCLBUF (NSYSTO, 4 ) -SCLBUF (N1ZR01, 4 ) )./ 
5 & (SCLBUF(NSYSTO*4)-29,5) )+0,5 

FEATUR (14+1 FEAT ) = FSYSTO/POSTDI 
C Ce qui suit reflate la presence ou 1" absence d'une 

"fenetre systolique" . 

APPELER EXTBYT ( SCLBUF (1,1), NPTS , NSYSTO- 12 ,20, IMIN , IMAX) 

10 FEATUR ( 1 5+IFEAT ) = SCLBUT(IMAX,l)-30 ! 3* dB inferieur. 

APPELER EXTBYT ( SCLBUF (1,3), NPTS , NSYSTO-12 ,20, IMIN , IMAX ) 

FEATUR ( 16+IFEAT ) = SCLBUF (IMAX, 3 ) -30 ! 9 dB inferieur. 

C Caract£ristiques bashes sur l'elargissement spectral 

(il y en a actuellement 28). 
FAIRE 750,1=1,4 t 71 
15 750 FEATUR ( IBR01+I+0+IFEAT ) = 

& FADIFB(LOCBUF,SCLBUF(1,I),NPTS,20,5) 

FAIRE 900, 1=1,4 ! 71 

900 FEATUR ( IBROl+I+4+IFEAT ) = 

& FADIFB (LOCBUF, SCLBUF (1,1), NPTS, NSYSTO, 5) 

20 FAIRE 1000, 1=1,4 1 75 

1000 FEATUR ( IBROl+8+I+IFEAT ) = 

& FADIFB { LOCBUF , SCLBUF (1,1), NPTS , NS YSTO+2 0 , 5 ) 

FAIRE 1600, 1=1,4 ! 79 

1600 FEATUR (IBR01+12+I+IFEAT)= 

25 & FADIFB ( LOCBUF , SCLBUF (1,1), NPTS , NSYSTO+40 , 5 ) . 

FAIRE 1700, 1=4 ! 83 

1700 FEATUR (IBR01+16+I+IFEAT) = 

& FADIFB(LOCBUF,SCLBUF(1,I),NPTS,120,15) ! 0,2 sec. 

apr&s ECG. 

30 C 87:0, 2 sec. aprfes ECG, largeur 3 

et 9 dB totales. 

FEATUR ( IBROl+21+IFEAT ) = 

& FADIFB(SCLBUF(1,1),SCLBUF(1,2),NPTS,120,20) 

» largeur 3 dB 

3 5 FEATUR ( IBROl+2 2 +IFEAT) = 

& FADIFB(SCLBUF(1,3),SCLBUF(1,4),NPTS,120,20) 

! largeur 9 dB 

C 89 



52 



2591884 



FEATUR ( I BR01 + 2 3+IFEAT ) = 

& FAD I FB ( SCLBUF (1,1), SCLBUF (1,2) ,NPTS ,10,10) 

I largeur 3 dB 

FEATUR ( I BROl + 2 4+1 FEAT ) = 
5 & FADIFB(SCLBUF(1,3),SCLBUF(1,4),NPTS,10,10) 

I largeur 9 dB 

FAIRE 1900, 1=1,4 ! 91 
1900 FEATUR ( I BROl + 2 4 +IFEAT ) = 

& FADIFB ( LOCBUF, SCLBUF (1,1), NPTS ,240,20) 
!0 ! 0,5 sec. apres ECG 

C DDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDDD 

C I FEAT: = nombre de caracteristiques actuelles, et est 
C egalement laiss€ en Evidence a la fin de la liste des 
C caracteristiques. 
15 5000 IFEAT= I FEAT+NUMFEA 

SI (NUMFEA .NE. (IBROl+28)) 

& IMPRIMER *, "les caracteristiques definies et engen- 

drees par BMAKI1-W ne correspondent pas". 

FIN 
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SOU S - PROGRAMME FEATUR ( LOCBUF , SCLBUF , NPTS , FEATS , NFEAT , 
SCALEF,NXO,NYO) 

COMMUN/FFTBUF/IREAL (1 024 ) , IMAG< 1024 ) 

COMPLEX CBUFF(l) 
5 ENTIBRM RMAG(l) 

EQUIVALENCE ( I REAL ( 1 ) , CBUFF { 1 ) r RMAG ( 1 ) ) 

OCTET LOCBUF { 1 ) , SCLBUF ( NPTS , 1 ) FEATS ( 1 ) 

OCTET SCR1(400),SCR2(400),SCR3(400) 
C Initialiser les coordonnees de tragage. 
10 NXLO= NXO+NPTS+10 

NYLO= NYO-30*SCALBF 
C Reproduire les donnees d 1 emplacement sur brouillon 1. 

LOCBUF (41)= LOCBUF (40) ! Eliminer le pic ECG 

APPELER STCOPY ( LOCBUF , SCR1 , NPTS ) 
15 C Calculer la premiere derivee et laisser dans le brouil- 
lon 1. 

APPELER DERIV(SCRl,NPTS,5,SCR2) 
APPELER TYPBUF(SCR2,NPTS) 
APPELER LISSE(SCR2,NPTS,4,SCR1) 
20 APPELER TYPBUF(SCR1,NPTS) 

APPELER SCABYT(SCR1,NPTS,2,5) 
APPELER ADDBYT ( SCR1 , NPTS ,50) 

APPELER PLTBYT ( SCR1 , NPTS , NXLO , NYLO , 1 • , SCALEF ) 
C Calculer les extremes dans la premiere moitie de la 
25 derivee. 

APPELER STCOPY ( SCR1 , SCR2 ,140) 

APPELER EXTREM ( SCR2 , 1 4 0 , 5 r MIN , MINI , MAX , MAXI ) 
C Calculer la FFT sur les donnees d* emplacement. 
C Reproduire les donnees d 1 emplacement dans la rangee 
30 d' entree FFT complexe. 

FAIRE 100, 1=1,512 
CBUFF(I)= 0 

SI (I .GT. NPTS) ALLER A 100 
CBUFF(I)= CMPLX(FLOAT(LOCBUF(I)- 30), 0.) 
35 100 CONTINUER 

C Obtenir le spectre de grandeur dans la rangee RMAG. 

APPELER MAGS PC ( NPTS , RMIN , RMAX ) 
C Avancer les coordonnees d f origine de la courbe. 
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NXLO= NXLO+NPTS+20 

Tracer la courbe du spectre de grandeur reel. 

APPELER SCAREA { RMAG , NPTS ,100, /RMAX ) 

APPELER PLTREA { RMAG , NPTS/2 , NXLO , NYLO , 2 . , SCALEF ) 

Tracer la courbe de l'intervalle de frequence pertinent 

uniquement. 

NYLO= NYLO- SCALEF*100-30 

APPELER PLTREA ( RMAG , NPTS/8 , NXLO , NYLO , 8 . , SCALEP ) 
Tracer la courbe logarithmique du spectre de grandeur. 
NYLO 53 NYLO- SCALBF*100- 30 
APPELER LOGRBA ( RMAG , NPTS ) 

APPELER SCAREA ( RMAG 7 NPTS r 5 0./ALOG( RMAX) ) 
APPELER ADDREA( RMAG, NPTS/2 ,40.) 

APPELER PLTREA ( RMAG , NPTS/2 , NXLO , NYLO , 2 . 9 SCALEF ) 

RETOUR 

FIN 
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REVEND I CATI ONS 

1. Procede non envahissant permettant de diagnostiquer 
automatiquement et object! vement le degr£ de st£nose d'un 
vaisseau a l'interieur du corps d'un patient , caract&rise 

5 en ce qu'il consiste : 

(a) a envoyer un signal ultrasonore convenable dans 
le corps , ce signal etant concentre au voisinage du centre 
du vaisseau pour mesurer 1 1 gcoulement du sang au milieu du 
vaisseau par la methode de variation de frequence due a I'effet 

10 Doppler ; 

(b) a recevoir le signal retrodiffuse de variation de 
frequence Doppler qui est engendre par I'ecoulement a l'int6- 
rieur du vaisseau ; 

(c) a enregistrer simultan£ment 1 ' electrocardiogramme 

15 (ECG) du patient avec le signal Doppler pour fournir un signal 
de minutage periodique permettant de separer le signal Doppler- 
en segments periodiques ; 

(d) a synchroniser le signal d'ECG avec le signal Doppler 
pour fournir un signal Doppler periodique normal i si de segments 

20 spectraux separ^s ; 

(e) a effectuer une analyse de frequence sur chaque 
segment spectral du signal Doppler par une methode de trans- 
formation de Fourier en temps reel, l 1 analyse etant d^clenchSe 
a partir d'un point predetermine du signal d 1 ECG r 1' analyse 

25 fournissant des donnees de frequence et d' amplitude ; 

(f ) a traiter les donnees de frequence et d' amplitude 
pour comprimer ces donnees, la compression comprenant l'obten- 
tion d'une moyenne d' ensemble des segments spectraux ; 

(g) a comparer la moyenne d' ensemble comprimee a une 
30 base de donnees de spectres connus et classes a I'avance 

contenus dans sa memoire, pour classer la moyenne d f ensemble 

en l'une de plusieurs categories indiquant le degre de stenose ; 

et 

(h) a diagnostiquer le degre de stenose du patient a" 
35 partir de la cat^gorie de classification. 

2. Procede selon la revendication 1, caracterise en 

ce que la comparaison de la moyenne d' ensemble comporte une 
reconnaissance de formes entre les donnees du patient et 
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les spectres de la base de donnees pour etablir la categorie 
de la moyenne d' ensemble en fonction des spectres connus 
contenus dans la base de donnees. 

3. Procede selon la revendication 2, caracterise en 

5 ce que la comparaison de la moyenne d 1 ensemble comporte une 
s£rie de decisions dans un format de decision hierarchique . 

4. Procede selon la revendication 3, .caracterise en 

ce que le format de decision hierarchique comporte au moins 
les decisions separees suivantes : 
10 (a) ma lade ou normal ; 

(b) si normal, stenose superieure ou inf^rieure a 50 %. 

5. Procede selon la revendication 1, caracterise en 

ce que le traitement des donnees de frequence et d" amplitude 
comprend les operations suivantes : 
15 (a) elimination des arythmies des segments spectraux 

analyses ; et 

(b) raise a la moyenne des segments spectraux d* environ 
15 k 20 segments par ailleurs successifs, sans arythmies, 
pour developper la moyenne d 1 ensemble. 

20 6. Procede selon la revendication 5, caracterisd en 

ce que le traitement comprend une operation de compression 
des donnees de la moyenne d* ensemble par selection d'un nombre 
pr^dSterminS de contours d 1 amplitude en fonction du temps, 
h des frequences predeterminees, dans la moyenne d 1 ensemble, 

25 en temps que representation de toutes les donnees de la moyenne 
d* ensemble. 

7. Procede selon la revendication 6, caracterise en 

ce que le traitement comprend le contour d' amplitude en fonc- 
tion du temps pour la frequence d f amplitude maximale et au 
30 moins un contour d" amplitude en fonction du temps de frequence 
superieure et de frequence inferieure pour representer l'6tale- 
ment du signal de variation de frequence Doppler veritable. 

8. Procede selon la revendication 7, caracterise en 

ce qu'on n' utilise pas plus de cinq contours pour la repre- 
35 sentation, ces cinq contours etant essentiellement le contour 
d' amplitude en fonction du temps pour la frequence d* amplitude 
maximale et deux contours d'echelle predetermines de part 
et d'autre du contour de frequence d'amplitude maximale. 
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9. Procede selon la revendication 5, caracterise en 
ce qu'on elimine les arythmies sur la base du fait que la 
durSe du signal d'BCG pour une arythmie est plus grande que 
la duree d'un segment type, presume normal, d'un temps supe- 

5 rieur a un temps predetermine. 

10. Proced6 selon la revendication 5, caracterise en 
ce qu'on eiimine les arythmies si la dur6e d^passe un temps 
predetermine, ou bien si la dur^e d^passe d'un temps predeter- 
mine la dur6e immediatement precedente. 

10 11. Procede selon la revendication 1, caract6rise en 

ce que le traitement comprend une operation de compression 
des donnees de la moyenne d' ensemble par selection d'un nombre 
predetermine de caracteristiques spectrales dans la moyenne 
d' ensemble en tant que representation de toute la moyenne 

15 d' ensemble. 

12. Procede selon la revendication 1, caracterise en 
ce que le vaisseau est une artfere carotide. 

13. Procede selon la revendication 1, caracterise en 

ce que la comparaison de la moyenne d' ensemble avec la base 

20 de donnees comporte une validation croisee de la classification 
pour chaque decision. 

14. Procede selon la revendication 13, caracterise en 
ce que la comparaison fait appel 4 une methode de resolution 
hyperplane de la reconnaissance des formes bas6e sur une 

25 analyse de regression lineaire, * la manifere de Pietrantonio 
et Jurs, communement appeiee Pjurs, selon laquelle I'algo- 
rithme de reconnaissance des formes minimise par iteration 
la fonction : 

K 

30 |Y - tanh (W Q + I wiXi)] 

i=l 

ou : Y = +1, si l'on sait que les donnees sont dans la cy=ite- 
gorie 2 

Y = -I, si l'on sait que les donnees sont dans la cate- 
35 gorie 1 

Y = la categorie predite 

K = nombre de caracteristiques 
wi = poids de la caracteristique 
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Xj[ = valeur mesuree de la caracteristique i 
w D = constante. 

15. Procede selon la revendication 1, caracteris6 en 
ce que le signal est un signal ultrasonore combing en mode 

5 B et a impulsions, fonctionnant a 5 MHz environ. 

16. Procede selon la revendication 1, caracterise en 

ce que 1' analyse de frequence est effectu^e dans un analyseur 
de spectres automatique, et le traitement et la comparaison 
des donnees sont effectu£s dans un ordinateur num£rique conve- 
10 nablement programme. 

17. Proc£d6 selon la revendication 1, caract6ris£ en 
ce que la comparaison comporte une reconnaissance de formes 

par rapport a une base de donnees enregis trees dans une meraoire, 
comprimees de fagon similaire avec les donnees de la moyenne 
15 d' ensemble, et non par rapport aux spectres complets. 

18. Procdd6 non envahissant permettant de diagnostiquer 
automatiquement et objectivement le degre de st^nose des 
art&res carotides d'un patient, caracterise en ce qu'il con- 
siste : 

20 (a) a envoyer un signal ultrasonore Doppler combine 

en mode B et a impulsions, fonctionnant a 5 MHz environ, 
dans le corps du patient, au voisinage de l'art&re carotide 
a examiner, ce signal etant concentre au voisinage de l'axe 
de l'art&re pour mesurer I'dcoulement central du sang a I'inte- 

25 rieur de I'artfere par la methode de variation de frequence 
due a I'effet Doppler ; 

(b) a recevoir le signal Doppler retrodiffusS qui est 
engendr£ par I'^coulement a l'int^rieur de l'artfere ; 

(c) a enregistrer simultan&nent l'ECG du patient pour 

30 fournir un signal de minutage periodique permettant de s£parer 
le signal Doppler en segments spectraux periodique s ; 

(d) k numeriser le signal Doppler dans un analyseur 
de spectres en temps reel ; 

(e) a definir les segments spectraux du signal Doppler 
35. dans 1' analyseur a partir d'un declenchement sur l'onde R 

de l'ECG, chaque segment spectral ayant une etendue suffi- 
sante pour permettre la detection d'un ecoulement direct, 
inverse ou instable a I'interieur de I'art&re ; 
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(f) a analyser chaque segment spectral dans l'analyseur 
par une methode de transformation de Fourier pour en extraire . . 
les donnees de frequence, d'amplitude et/ou de temps ; 

(g) a transferer les donnees de frequence, d'amplitude 
et/ou de temps de l'etape (f ) dans un ordinateur num^rique 
convenablement programme , ainsi que les donnees du signal 
d ' ECG ; 

(h) a calculer une moyenne d 1 ensemble des segments spec- 
traux en filiminant les arythmies de la moyenne, sur la base 
de la plus grande dur£e des arythmies par rapport aux seg- 
ments types, presumes normaux, et a faire la moyenne d' envi- 
ron 15 a 20 segments spectraux successifs, par ailleurs accep- 
tables, a l'aide de la formule : 

K 

15 E (f,t) I E (f,t,n)/K 

n=l 

dans laquelle : 

E est I 1 amplitude de la transformation de Fourier a 
f , t et n ; 

K est le nombre total de segments dans la mpyenne ; 
f est la frequence ; 

t est le temps (mesur6 par rapport a l'onde R de 1'ECG) ; 

et 

n est le numero d'ordre du segment ; 

(i) a comprimer la moyenne d' ensemble pour obtenir une 
representation spectrale convenable pour la classification 
par reconnaissance de formes en (1) choisissant la frequence 
d f amplitude maximale de la moyenne d' ensemble en temps qu' indi- 
cation du signal veritable en l'absence de bruit, (2) crSant 
un contour d'amplitude en fonction du temps a cette frequence 
pour la duree du segment et pour cette frequence, et (3) 
calculant un £talement representatif du signal . en incluant 
un ou plusieurs contours d'amplitude en fonction du temps 
au-dessus et au-dessous du contour de frequence d'amplitude 
maximale, la moyenne d' ensemble comprim^e representant au 
moins trois contours d'amplitude' en fonction du temps ; 

( j ) a comparer la moyenne d 1 ensemble comprim^e a une 
base de donnees dans la memoire de laquelle se trouvent des 
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donnees semblables de contours de moyennes d' ensembles prove- 
nant de vaisseaux px*esentant des degres connus de stenose, 
la base de donnees etant acquise a partir de donnees de varia- 
tion de frequence Doppler sur des patients dont le degre* 

5 de stenose a ete confirme* par arteriographie de contraste, 
la comparaison s'effectuant par reconnaissance de formes 
dans un processus de. decision hierarchique utilisant au 
moins les e tapes de decisions suivantes : 
(1) vaisseau bouche ou ouvert ; 

10 (2) vaisseau malade ou normal, s'il est ouvert ; et 

(3) stenose sup£rieure ou inf£rieure a 50 %, si le vais- 
seau est malade i 

cette comparaison permettant de classer la moyenrie d f ensemble 
dans une cat6gorie k l'aide d'une methode de resolution hyper- 
15 plane ; et 

(k) a afficher la classification en tant que diagnostic 
du degre de stenose. 

19. Procede selon la revendication 18, caracterise en 

ce que le signal Doppler provient d'au moins I'artere carotide 
20 primitive basse et de I'artfere carotide interne proche. 

20. Procede selon la revendication 19 $ caracterise en 
ce que la methode de resolution hyperplane est basee sur 
une analyse de regression lin^aire it la maniere de Pietr- 
antonio et Jurs, communement appelee Pjurs, selon laquelle 

25 l'algorithme de reconnaissance des formes minimise par itera- 
tion la fonction : 

K 

I Y - tanh (W 0 + I WjXj.) J 
i=l 

30 ou : Y = +1, si l'on sait que les donnees sont dans la cat£- 
gorie 2 

Y » -1, si l'on sait que les donnees sont dans la cate- 

gorie 1 

Y est la categorie predite 

35 K est le nombre de caracteristiques 

W£ est le poids de la caracteristique 
Xi est la valeur mesuree de la caracteristique i 
w Q est une constante. 
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21. Proc£de selon la revendication 20, caracterise* en 
ce que la classification dans une cat£gorie a chague niveau 
de decision hi^rarchique est confirmed par une analyse de 
validation croisee. 

22. Procede* selon la revendication 21, caract6ris£ en 
ce que le proc6d<5 de validation croisee est un precede dit 
"demi-echantillon plus complement". 

23. Proc6de* selon la revendication 22, caracteris£ en 
ce que les caract^ristiques h comparer entre les donndes 
de variations de frequence Doppler et la base de donnees 
comprennent les caracteristiques suivantes pour chaque deci- 
sion : 

(a) vaisseau malade ou normal ; 

(b) si le vaisseau est malade, stenose sup^rieure ou 
15 infe*rieure & 50 % ; 

(c) si la stenose est infSrieure & 50 % f stenose comprise 
entre 50 et 35 % ou inf^rieure & 35 % ; et 

(d) si la stenose est infSrieure a 35 %, stenose comprise 
entre 35 et 15 % ou inf^rieure & 15 %. 

20 24 • Proceae* selon la revendication 18, caract^rise en 

ce que la decision entre vaisseau occlus ou ouvert comporte 
des essais pour determiner s'il y a un signal de variation 
de frequence Doppler en provenance de 1'artere carotide interne 
et de 1'artere carotide primitive pour verifier que les paro- 
tides sont ouvertes et non completement bouchees. 

25. Procede non envahissant permettant de diagnostiquer 
automatiquement et objectivement le degre* de stenose des 
arteres carotides d'un patient humain h partir des donnees 
de variation de frequence Doppler & ultrasons en provenance 
de 1'artere carotide primitive basse et de 1'artere carotide 
interne proche, caracterise en ce qu'il consiste : 

(a) a envoyer un signal ultrasonore Doppler combine* 
en mode Beta impulsions , fonctionnant & 5 MHz, dans le 
corps du patient, au voisinage de 1'artere carotide h exa- 
miner aux endroits voulus, ce signal 6tant concentre* au voisi- 
nage de I'axe de 1'artere pour mesurer l'Scoulement central 
du sang a l'interieur de 1'artere en chaque endroit par le 
proced^ de variation de frequence Doppler ; 



25 
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(b) a recevoir le signal Doppler rStrodiffuse qui est 
engendr£ par l'ecoulement a I'interieur de l'artere ; 

(c) a enregistrer simultan&nent l'electrocardiograrame 
du patient pour fournir un signal de minutage periodigue 

5 permettant de separer les signaux Doppler en segments p6riq- 
digues correspondants a des cycles cardiaques successifs ; 

(d) a numeriser le signal Doppler dans un analyseur 
de spectres a transformation rapide de Fourier discrete en 
temps r£el ; 

10 (e) a d6finir des segments spectraux du signal Doppler 

dans 1' analyseur, sur la base d'un declenchement sur l^nde 
R de I'ECG ; 

(f ) a analyser chague segment spectral dans 1" analyseur 
par une methode de transformation de Fourier pour en extraire 

15 des donnees de freguence, d' amplitude et/ou de temps ; 

(g) a transferer 1' analyse de l'6tape (f ) dans un ordi- 
nateur numerique convenablement programme/ ainsi gue les 
donnees du signal d'ECG ; 

(h) a calculer une moyenne d' ensemble des segments spec- 
20 traux en chaque emplacement en elirainant les arythmies de 

la moyenne sur la base d'une dur£e plus longue des arythmies, 
et en faisant la moyenne d 1 environ 20 segments spectraux 
successifs, par ailleurs acceptables, a l'aide de la formule : 
K 

25 E (f,t) = I E (f,t,n)/K 

n=l 

dans laguelle : 

E est 1' amplitude de la transformation de Fourier a 
f , t et n ; 

30 K est le nombre total de segments dans la moyenne ; 

f est la frequence ; 

t est le temps (mesurd par rapport a I'onde R de l'ECG) ; 

et 

n est le num^ro d'ordre du segment ; 
35 (i) a comprimer la moyenne d' ensemble de I'ecoulement 

en chague emplacement pour obtenir une representation spec- 
trale convenable pour une classification ef f icace et simplifiee 
par reconnaissance de formes, en (1) trouvant la frequence 
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d' amplitude maxima le de la moyenne d' ensemble a partir d'une 
courbe d 1 amplitude en fonction de la frequence pour la moyenne 
d 1 ensemble du cycle cardiague complete la frequence d' amplitude 
maximale etant une indication du signal veritable en 1» absence 

5 de bruit, (2) errant un contour d f amplitude en fonction 
du temps a cette frequence pour la totalite du cycle cardiaque 
du segment, et (3) calculant un £ta lenient representatif 
du signal en incluant les contours d" amplitude en fonction 
du temps de 3 dB et 9dB au-dessus et au-dessous du contour 

10 de frequence d' amplitude maximale, les donnees comprimees 
du contour de moyenne d 1 ensemble etant essentiellement cinq 
contours d' amplitude en fonction du temps ; 

(j) a comparer des caracteristiques sel'ectionn£es de 
la moyenne d v ensemble comprimee £ une base de donnees conte- 

15 nant dans sa meraoire des caracteristiques semblables selec- 
tionnees parmi des donnees de contour de moyenne d' ensemble 
obtenues dans des conditions connues de diagnostic prealable, 
la base de donnees etant acquise a partir de donnees de varia-? 
tion de frequence Doppler sur des patients dont le degre 

20 de stenose a et6 confirm^ par art^riographie de contraste, 
la comparaison s'effectuant par reconnaissance de formes 
dans un processus de decision hierarchique utilisant au 
moins les etapes de decision suivantes : 
(1) vaisseau bouch£ ou ouvert ; 
25 (2) vaisseau ma lade ou normal, s'il est ouvert ; et 

(3) stenose sup£rieure ou inferieure h 50 %, si le vaisseau 
est malade, 

la comparaison permettant de classer la moyenne d' ensemble 
du patient dans une categorie par 1 ■• interm^diaire d'une m£thode 
3 Ode resolution hyperplane qui est basee sur une analyse de 
regression lin^aire k la manidre de Pietrantonio et Jurs, 
communement appelde Pjurs, selon laquelle -J:* algorithms de 
reconnaissance des formes minimise par iteration la fonction : 

K 

35 [Y - tanh (W Q + £ w ± Xi) J 

i=l 

ou : Y = +1, si Ton sait que les donnees sont dans, la cate- 
gorie 2 
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Y s -1, si l'on salt que les donnees sont dans la cat£- 

gorie 1 

Y est la categorie pr£dite 

K est le nombre de caracteristiques 
5 w£ est le poids de la caract£ristique 

Xi est la valeur mesuree de la caracteristique i 
W 0 est une constante 
la methode utilisant jusqu'a cinq caracteristiques pour chaque. 
decision ; 

10 (k) & operer une validation crois^e de chaque decision 

a I'aide d'une methode iterative "derai-echantillon plus comple- 
ment" pour verifier la decision et pour reduire la variance 
inherente aux donnees ; et 

(1) k afficher la classification determine par 1» analyse 

15 par ordinateur en tant que diagnostic du degr£ de stenose. 

26. Proc^de non envahissant permettant de diagnostiquer 
automatiquement et objectivement le degre de stenose des 
arteres carotides d'un patient humain & partir de donnees 
de variation de frequence Doppler ultrasons en provenance 

20 de l'artere carotide primitive basse et de l'artfcre carotide 
interne proche, caracterise en ce gu'il consiste : 

(a) h produire un signal de variation de frequence 
Doppler h chaque emplacement, incluant environ 20 h 25 cycles 
cardiaques successifs ; 

25 (b) h enregistrer simultan&nent 1 •Slectrocardiogramme 

du patient pour f ournir un signal de minutage permettant , • 
de separer les signaux de variation de frequence Doppler 
en segments spectraux periodiques correspondants a des cycles 
cardiaques successifs ; 

30 ( C ) a faire la moyenne des signaux provenant des segments 

de fagon h produire un signal de variation de frequence Doppler 
repr£sentatif de I'ecoulement £ I'interieur de I'artfere, 
h l 1 emplacement correspondant ; 

(d) a comparer la moyenne d* ensemble de l'etape (c) 

35 avec une base de donnees de moyennes d* ensemble classees 

h l'avance, correspondant a une stenose connue, sur un ordi- 
nateur numerique, par une methode de reconnaissance de formes 
pour classer le degre de stenose ; et 
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(e) a afficher la classification determine par la compa- 
raison par ordinateur en tant que diagnostic du degre de 
stenose. 

27. Dispositif permettant de diagnostiquer automatic 
5 quement et objectivement le degr6 de stenose d'un vaisseau 

& partir de donnees Doppler combines en mode Beta impulsions, 
analysees en segments spectraux dans un analyseur de spectres 
a transformation rapide de Fourier en temps r6el, caract6ris£ 
en ce qu'il comprend un ordinateur num£rique convenablement 
10 programme pour effectuer la reconnaissance de formes des 
segments spectraux ou d'une moyenne d 1 ensemble des segments 
spectraux pour classer le degx6 de stenose en fonction d'une 
base de donnees de stenoses connues en provenance de segments 
spectraux classes k I'avance. 
15 28. Dispositif permettant de diagnostiquer automatique- 

ment et objectivement le degr<§ de stenose d'un vaisseau humain 
par une technique non envahissante, caracterise en ce qu'il 
comprend : 

(a) un scanographe Doppler a ultrasons a fonctionnement 
20 combing en mode B et en impulsions, pour envoyer un signal 

concentre au centre du vaisseau et pour recevoir le signal 
de variation de frequence Doppler r6trodiffus6 repr£sentatif 
de l'ecoulement sanguin a l'int&rieur du vaisseau ; 

(b) un analyseur de spectres pour transformer le signal 
25 de variation de frequence Doppler analogique en un signal 

numerique frequence-amplitude-temps ; 

(c) un detecteur d'onde R d'dlectrocardiogramme pour 
minuter simultanement le signal de variation de frequence 
Doppler en segments p^riodiques correspondants & des cycles 

30 cardiaques successifs ; 

(d) un ordinateur num&rique, comportant une base de 
donnees dont la m&noire contient des segments, classes a 
l'avance, de signaux de variation de frequence Doppler corres- 
pondants k une stenose connue, permettant de comparer le 

35 signal du patient k la base de donnSes dans une methode de 
reconnaissance de formes pour classer le degre de stenose ; 
et 

(e) un visuel pour indiquer la classification d^termin^e 
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dans I'ordinateur numerique. 

29. Dispositif selon la revendication 28 , caracterise 
en ce que l'analyseur de spectres fournit a l'ordinateur 
une transformation rapide de Fourier discrete et numerisee 

5 des donnees de variation de frequence Doppler. 

30. Dispositif selon la revendication 29, caracterise 
en ce que l'ordinateur compare le signal du patient h la 
base de donnees en : 

(a) eliminant les ary thirties du signal ; 
10 (b) produisant une moyenne d 1 ensemble de segments perio- 

diques successif s du signal identif i6 par une position de 
minutage sur l'onde R. dans le signal d'ECG ; 

(c) comprimant le signal & au moins un contour d f ampli- 
tude en f onction du temps repr^sentatif du veritable ecoule- 

15 ment sanguin en 1' absence de bruit ; et 

(d) rangeant le degrd de st£nose du patient dans une 
categorie en comparant les caract£risques s£lectionnees du 
signal comprint k des caract^ristiques semblables extraites 
de la base de donn£es de signaux provenant de vaisseaux 

20 atteints de stenoses connues. 

31. Proc6d£ permettant de diagnostiquer le degr6 de 
stfinose d'une artfere carotide par une technique ultrasonore 
non envahissante, caract£ris£ en ce qu'il consiste a comparer 
les caract£ristiques extraites des donnees de variation de 

25 frequence Doppler de plusieurs battements cardiaques successif s 
du patient & une base de donnees de caract£ristiques Doppler 
analogues en provenance de vaisseaux atteints d'une stenose 
connue, selon un format de reconnaissance de formes, pour 
identifier le degr£ de st&iose. 

30 32. Proc6d£ selon la revendication 31 r caract£ris6 en 

ce que le diagnostic comprend la prise de decisions catego- 
rielles du degr£ de stenose dans un format hierarchique. 

33. Proced£ selon la revendication 32, caracterise en 
ce qu'on fait la moyenne des donnees des battements cardia- 

35 ques successifs pour comprimer les donnees en une moyenne 
d 1 ensemble . 

34. Dispositif permettant de diagnostiquer automati- 
quement et object! vement le degre de stenose d'une artere 
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carotide d'un patient par 1 * intermediate d'une technique 
ultrasonore non envahissante, caract£ris6 en ce qu'il com- 
prend : 

(a) un scanographe-recepteur a ultrasons ; 
5 (b) un moyen permettant de comparer les caract£ristiques 

extraites des donnees ultrasonores de plusieurs battements 
cardiaques successifs du patient a une base de donnees de 
battements cardiaques analogues provenant de vaisseaux atteints 
.d'une stdnose connue, de fagon k placer les donnees du patient. 
10 dans une categorie de stenose ; et 

(c) un visuel servant k presenter la categorie du degr6 
de stenose determin£e a partir des donnees ultrasonores. 
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